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ABSIRACT 



Expeiisents were performed on laser produced plasma 
from aluminum target discs in a vacuum chamber at 



-6 

pressures a few times 10 torr. Plasma analysis was 

achieved using a floating double electrostatic probe 
of tungsten wires biased at -15 volts DC. In 
analyzing 376 oscilloscope trace pictures of the probe 
signal/ four velocity signals occurred with 
regularitj. Ihe plasma velocities detected by the 
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probe were determined to be 3.7X10 cm/sec, 1-5X10 
6 

cm/S€C/ 5.6X10 cm/sec (the bulk of the plasma), and 

e 8 

14.0X1C cm/sec. A fifth, very fast signal at 1.0X10 

cm/sec was also detected but it was not consistently 

present. Investigation of the velocity distributions 

as a function of time showed them not to be Boltzmann 

distributions, tut exponential power distributions . 

Investigation cf the crater formed by each laser 

pulse resulted in the determination of the following 

crater dinersicns; crater depth 0.01 cm/shct, crater 

diameter cf 0.11 cm and a total spot size diameter of 

0.70 cm for the laser-target interaction. 
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I. 



INTflOEDCTIOM 



laser plasma production has been of interest to 
researchers since the early days of laser study. The use of 
the laser as a scientific instrument quickly dre% the 
attention of the academic community. Laser-target 
interaction was soon a topic of considerable research 
effort. Industrial and military applications of 
laser-dominated technology are considerable. Most recently 
it has beer employed in research to develop the fusion 
reactor as an efficient, controllable source cf energy, tc 
develop a rapid, secure communication system, and to develop 
the much publicized laser weapon systems of future warfare. 
The effects cf high power laser interaction with matter 
include <1) beatinc £42], (2) melting, (3) vaporization of 
solid materials, (4) emission cf charged particles, (5) 
emission cf neutral particles [54,11], (6) plasma 
production, (7) electrical discharge in gases produced by 
high energy laser interaction, and (8) application cf these 
effects tc material processing. 

Addressing laser interaction with solid materials, 
particularly metals, the mechanism cf the interaction has 
been the subject of research for over ten years. It is now 
uncerstocd that in metal targets struck by a laser beam, 
some of the light energy is absorbed by the electrons and a 
large part is reflected, the amount of reflection being 
dependent cn the target material [42]. The amcurt of 
radiation energy and the time-length cf the laser pulse 
striking the target are the controlling parameters that 
.determine the results of the interaction. 

The electrons in the conduction band will be excited 
to a higher energy level as incident photons are absorbed. 
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Tb€ energy ficm these interactions is transferred through 
ccllisiocs in the target material to other electrons and 
ions. The mean collision time for the electrons is 
one-tenth picosecond, so it is apparent that many collisions 
can occur in the time of a laser pulse measured in 
nanoseconcs. The photon energy is, for all practical 
consideraticES , instantly turned to heat through particle 
collisions, rapidly heating a small quantity of material to 
very high temperatures. A laser-target interaction of 
sufficient energy will produce a plasma that is a dense 
cloud of ionized particles. Plasma production occurs at 

9 2 

laser irradiances of the order of 10 W/cm and greater 
[42]. 

The early response of the laser produced plasma has 
marked characteristics that are in need of further study. 
The plasma plume develops guickly after the laser pulse 
strikes the target. The plume development is at least 
partially described by computer analysis at lavrence 
Livermore laboratories [32], Sandia Laboratories, and Naval 
Research laboratory [37,38]. The Naval Research Laboratory 
in fiashington, E.C. has a computer model based on a coupled 
system of hydromagnetic equations in cylindrical coordinates 
including self-generated field effects, rate equations, 
radiation transport, frequency diffusion of incident photons 
and radiation pressure [37]. 

The early response of the laser produced plasma which 
precedes the formation of the main plume is more difficult 
to explain [52]. It is hypothesized through experimental 
observation that there are four possible sources of fast 
propagating early signals: laser induced electron emission 
[52,62], fast plasma blowoff [24], electromagnetic signals, 
and ionization of the background gas by plasma produced 
photons [27] or ionizing potential waves [20]. 
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Ihis study was undertaken to develop further understanding 
of the origin of the early plasma emission through 
experimental study of the velocities of plasma features. 
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II. HECENl BELATED EESiAECH 



The E€WE€ss of laser produced plasma research has not 
jet produced manj comprehensive textbooks cn the subject. 
Most of the literature is in the scientific journals. The 
most comprehensive book, judging by the number of times it 
is referenced by other authors and the extensive 
bibliography it contains, is ref. [42] by J. F. Beady, It 
is five years old and provides reference to most of the 
published articles on laser effects prior to 1971. It 
contains two chapters directly related tc laser target 
interaction with a metal surface and several more chapters 
covering other target materials, i.e., gases, transparent 
material, semi-conductors, and non-metallic solids. Chapter 
Three, "Effects Caused by Absorption of Laser Radiation at 
Opague Surfaces," discusses investigations into heating 
without tber mcdyna aic phase change, laser induced melting 
and laser produced vaporization. 

Laser induced heating without phase changes is 
basically a discussion of temperature change using classical 
thermody cam ics formalism. The laser photons interact with 
the electrcrs in the solid and distribute their energy 
through the metal. The mean free time fcr collisions 
between electrons in a conducting metal is cf the crcer of 

- 13 

10 seconds. The laser radiation energy cf a Q-switched 

- 8 

or normal pulsed laser is slowly deposited (t=10 seconds 

as compared tc the electron mean free collision time of 

.13 

•10 seconds). A condition of thermal eguilibrium can be 

assumed tc exist. 

J. H. Eechtel also has investigated laser-induced 
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t€.ferat«ie changes of the absocbing eedia [8] 

Incident iitadiancee such that neither halting nor 

vatcriraticn cccurs. Here the heat-conducticn eguatior 



V J ( r , t ) + pc 3T (r , t ) - A ( r , t ) 

3 1 



is eolved and is called the volume generation model, in 
vhich I in a temferature function, j' is the thermal energy 
crossing a unit area per unit time, fc. is the heat capacity 
net unit volume and 1 is the net energy per unit volume per 
unit time generated eithin the solid. An alternative method 
of formulating the problem of laser heating of ahscrting 
material is called the surface generation model. In this 
solution the heat is assumed to be generated on the surface 
of the target and the homogeneous heat-conduction eguation 
is solved «ith loundary conditions at the irradiated 

surface. Eiperimental evidence is derived for CO^ 

radiation impinging on a tungsten target to agree uith the 

theory developed [8]. The surface generation model can be 
u=ed to predict laser-induced temperature changes produced 
b, pulsed lasers. Por very short pulsed, picosecond, 
mode-locked solid state and dye lasers the volume generation 
mcdel is acre accurate. 

The nc-phase-change laser heating problem is addressed 
by J. I. schrieicpf [39]. He solves the beat flo. eguataon 
for a seni-irfinite solid and a slab of finxte thickneo^ 
(thin film) and compares the results to experiaenta 

observations on aluminum and stainless steel. -*• 
and S. £. laff [21] have reported on laser heatxng of thxn 
films. Ihe iicblem is solved using the one dimensional heat 
■flow equation, a numerical example for a silicone fila cn a 
silica substrate illuminated with a CW argon-ion laser is 
compared vith experimental observations by Ghez and lafr. 
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Beady [42] continues in Chapter Three 

laser-induced melting and a theoretical solution usi 
heat equation. Melting or reducing a solid to its 
phase without vaporizing it requires a very fine 
division when using a laser energy source. The 
lasers ate not especially suited for that applicatio 
CH lasers which can be adjusted to variable energie 
well suited for melting. Schriempf [39] work 
seai-inf icite slab melting problem theoretically by u 
the sinple heat flow equation 

3^T ^ 3T 

9z^ ~ K 3t “ ° 



assuming the melted material is immediately removed, 
is a realistic assumption demonstrated by J. A. Fcx 
fcx showed that two lasers (one pulsed and one CH 
provide a workable technique tc clear the molten 
material awaj rapidly enough tc consider it as being r 
ianediately. Ghez and Laff [21] address laser meltin 
film and substrate as a three-phase Stefan problem of 
film, solid film, and substrate. The Stefan prooleir 
in which a given material transforms from cne thermcu 
phase tc ancther with the emission or absorption o 
[12]. Assumptions are made to remove anal 
difficulties that make the problem otherwise insoluble 
scluticn tc the problem is found by the simple heat 
equation as in Schriempf's work. 

Considerable attention is given tc laser-i 
target vaporization and plasma production in the scie 
lit-erature. It is probably the most important aspe 
laser-target interaction when consideration is given 
ever more energetic laser pulses and the shorter 
■ widths cf state-^cf-the-art lasers. Vaporization and 
prcducticc are of direct interest to scientists atte 
tc design fusion reactors who need tc knew exactly w 
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bapfening at the fuel surface when the laser beams strike 
the deuteriuB-tritium fuel pellet [29]. Vapcrizaticn cf the 
target material is very readily attained with pulsed and 
Q-switchec lasers. Eeady [42] emphasizes the experiaental 
neasurement s cf material removal and seme metheds of 
calculating it. There is a great difference between the 
behavior cf c target surface struck by a millisecond laser 
pulse ard a target surface struck by a nanosecond or 
picosecond pulse. Intense, short laser pulses dc not 
produce much vaporized material removing only a small amount 
of target material frem the target surface with high flux 
density. longer lower-powered pulses frem a normal pulsed 
laser on tie ether hand produce deep, narrow holes in the 
target. ixperimental data on the amount of material 
vaporized hy a laser pulse is not very reproducible. It is 
dependent or exact duplication of experimental conditions 
and even the laser pulse is only reproducible to within ±20^ 
cf the rated energy from pulse to pulse on many high pewer 
lasers. Seady [42] presents several graphs cf laser crater 

a 

depth versus laser flux density (H/cm ) and laser energy 

output (Joules). The heat flew eguation is used to give a 

quick, simple estimate of the depth of vaporization. is the 
target material vaporizes, the thermodynamic phase changes 
that take place invalidate the heat equation for all but 
rough estimates. The thermal conductivity, heat capacity, 
and density are all temperature dependent. Ihe reflectivity 
is temperature and thermodynamic phase dependent. The 
laser-target crater depth is therefore not actually a simple 
function cf position and temperature as the heat equation 
indicates . 

Cne model cf the processes involved in vaporization by 
a C~sttitchec laser considers the pressure produced tj the 
prcQucticn cf a small amount of blowoff material early in 
the laser pulse. This recoil pressure raises the boiling 
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tenperature cf the target atove the noraal atnjospheric 
pressure vapcrizaticn temperature. If then, one considers 
that the increase in vaporization temperature of the surface 
material is sufficiently high, the surface will be prevented 
ficn: vapctizing further and the material will continue to 
heat to a temperature above normal vaporization temperature 
as laser light continues to be absorbed by the target 
surface. Ike target surface will eventually reach the 
critical pcict where vaporization can occur. This mcdel 
gives reascrafcle agreement with experimental observations 

9 2 

for one range cf laser flux density (10 M/cm ). It makes no 

allcwances fcr shielding or partial shielding of the target 
surface ficm the laser beam at the target spot as the 
blowoff material (plasma plume) becomes hct, ionized, and 
opague. 

at sufficiently high laser power density the blowoff 
can be assumed tc be a plasma. Gas dynamics, plasma 
absorption mechanisms, thermal ionization processes and 
reradiaticn by the plasma are all included in a model given 
by Beady [^2j. This model predicts experimental facts with 
relatively good accuracy. Graphs for vapcrizaticn depth 

2 

(crater depth) as a function of energy density [J/cm ], 

velocity cf surface vaporization [cm/sec] vs. laser flux 
z 

density [fc/ca ], target surface temperature as a function of 
laser flux decsity and vaporization depth versus laser pulse 
duration are presented for comparisons. 

These ncdels each emphasize a particular aspect cf the 
interaction. They each give a reasonable treatment cf one 
of several pheromera, but there are no comprehensive models 
that compile all of the phenomena into a single solvable 
model . 

Continuing, Schriempf [39] presents a solution tc the 
problem cf uriform continuous laser irradiation of a sclid 
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taiget; the cne dimensional problem. The melted vclcme of 
the target is assumed to be in place until it reaches 
vaporization temperature, when it is assumed to evaporate 
ccmpletelj. The pulsed laser situation for laser radiation 
is assumed. This is not an unreal situation as demonstrated 
by fox £1S]. One is usually interested only in power 
densities below the breakdown threshold of air because at 
higher power densities the energy cannot reach the target if 
it is in the ambient atmosphere. The breakdown levels are a 
function of laser radiation wavelength, spot size, laser 
pulse length, laser pulse shape and the contaminants in the 
atmosphere (i.e,, smog, fog, chemicals, smcke and humidity). 
It is for this reason that mcst experiments ccncerning 
laser-target interaction are done under at least partial 
vacuum conditions. High power density la-ser pulses are 
considered by Schriempf £39] in laser-target interactions 
assuming that the plasma blowoff does not shield the target 
spot from the laser radiation. In addition to the thermal 
change to the target, there is a strong pressure built up on 
the target surface due to reccil from the vapor blcwcff. 
This is where the spall and deformation of the target 
surface is suspected to ccme from. A crater of much 
Sfflccther geciretrical shape can be produced with a two laser 
system which removes the melted material quickly before it 
can vaporize and set up a shock front [19]. A theoretical 
calculaticn is done for the semi-infinite slat and numerical 
examples are carried cut for titanium and aluminum [39*. 

B. i. Clstad and D. H. dander [40] in a study of the 
vaporization cf iron by a laser pulse, have used the 
transient teat conduction equation tc solve the 
laser-interaction problem. Two assumptions are made ir this 
wcrk: (1) heat conduction occurs only in a direction rormal 
tc the target surface, and (2) because metals are cpague to 
visible light, the laser radiation is represented by a 
surface teat source. The temperature and time histories of 
the solid duiing a pulse may be compared by scluticn of the 
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transient heat equation. First, the surface temperature 
rises, meltiig and vaporizing the target. Neutral atoas or 
mclecules are emitted into the gas or vacuum arcund the 
solid. If the surface temperature exceeds the melting point 
of the sutstrate, rapid vaporization of the melt occurs and 
the escaping vapor exerts a recoil pressure on the liquid 
which causes it to flow outward along the surface. This 
movement cf liquid is observable as a crater at the point of 
impact cf the laser radiation. The recoil pressure into the 

9 

solid tarcet, when using laser beam energy densities cf 10 
(B/cm ) and greater, is enough to eject small globules of 

molten metal in a direction away from the target surface. 

P. S. E. Bei, D. J. Nelson, and E. B. Ball [59] 
present a simple model for evaporation in which the surface 
temperature is raised according to its heat capacity, the 
absorbed energy and the heat conductivity to the bulk. When 
the vapor pressure becomes comparable to that cf the 
surrcundirg atmosphere, surface boiling begins. Ir the 
absence cf superheating, the evaporated particles velocity 
would obey the three-dimensional Boltzmann distributioc. 

For plasma production. Beady [42] and Huddlestcne and 
Leonard [25] both have a good review cf the different 
techniques used in experimentation to investigate the 
characteristics of the plasma plume in high energy 
laser-target interactions. Photographic techniques, 
interferometric techniques, optical spectroscopic studies. 
X-ray emission, charge collection measurements with shielded 
probes, pressure pulses, and magnetic field considerations 
are all briefly described by Ready and in greater detail by 
Huddlestcne with much attention to extensive bibliographies 
to support each method. 

laser- target interaction can be characterized by two 
conditiocs at the impact area. The first is a low 
temperature region where the target material is barely 
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iccized and the second a high temperature region where the 
target material is in a plasma state [301. In the low 

temperature region the target material is in the process of 
changing thermodynamic states from solid to vapcr. The 
phase -transitions are characterized by a strong, djnamic 
transition ficm one state to another in a short period of 
time. In the high temperature region a hot plasma is 
produced that expands to form the plasma plume in a 
direction away from the impact area. The velocity cf the 
expansion and the composition cf the expanding plasma are 
areas of much concern to researchers. 

Laser produced plasma particle velocities are the 
basis of this thesis. The following is a summary of seme of 
the research dene on laser produced plasma particles tc date 
and the findings using various diagnostic techniques. 
Demtrcder and Jantz [17] found, using a tetrode screen 
detector system, that the velocity of hydrogen ions in a 

7 

laser induced plasma blowoff was UIIO cm/sec with a SCO KH 

7 

incident laser power, carbon ions had a velccity of 1.6X10 
cm/sec fer a 200 MF incident laser power, and copper ions 

7 

had a velccity of 10 cm/sec also for incident power equal 

to 200 MF. & Q-switched ruby laser was used in these 
experiments with a pulse width of 30 nsec. The incident 

It 2 

pewer densities ranged up to 5X10 H/cm . The detecter was 

placed 170 cm frem the target. 

Incue, Kawasumi, and Miyamota [26] found plasma 

particle velocities an order of magnitude smaller for carbon 

iens usirq ar ion energy analyzer. A ruby laser was used in 

these erperiments with a maximum power cf 60 MW and total 

energy 1.2 0. The target was a polyethylene (CH ) 

2 n 

filament. The only species they detected in the plasma was 
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2 + 

C . Two velccities were determined (parallel velocitj and 

6 

lateral velccity) for the carbon ions, 3.4X10 cm/sec and 
6 

1.5X10 ca/sec. 

Gregg and Ihomas [23] used a Q-switched ruby laser 

10 11 2 

with power densities from 1.5X10 to 5X10 W/cm . They 

experimented using lithium hydride (liH) , lithium, 
beryllium, carbon, sulfur, aluminum, zinc, and silver 
targets. Fcr each target material a graph was plotted of 

2 2 

the leg cf tte ion velocity squared [cm /sec ] versus the 

2 

leg of flux density [W/cm ]. The ior detector was a 
phctcmultiplier tube 40 cm frem the target. An analytical 



relationship was 


discovered for 


the ion velocities of each 


target material. 


It was found by 


graphical analysis that 




2 


X 




(V ) 


= AE and. 




avg 






2 


y 




(V ) 


= BE 




max 




given E egcal 


to the peak Intensity of the laser pulse in 


2 


10 


1 1 


«/cm from about 


1.5X10 to 5X10 


and A, X, Cf and y are 



given by the following table. 
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Table 1. 



Material 


X 


y 


A 


B 








9 


e 


liH 


0.459 


0.540 


2.12X10 

9 


5.05X10 

8 


li 


0.425 


0.533 


3. 99X10 

8 


4.93X10 

e 


Ee 


0.521 


0.565 


4.14X10 

8 


2.72X10 

e 


C 


0.561 


0.619 


1.20X10 

8 


5.48X10 

7 


5 


0.469 


0.627 


7.48X10 

6 


2.72X10 

4 


ai 


0.695 


0.859 


1 . 90X10 

1 0 


5.82X10 

8 


Zn 


0.330 


0.500 


1.06X10 

9 


3.44X10 

9 


iq 


0.339 


0.377 


6.24X10 


5.24X10 


let an aluainum 


ta rget 


their findings for 


ion velocities 




6 


7 






ranged frem 5X10 


to 1 


.7X10 cm/sec. 




Beady £43] 


used 


visual micrcscopic analysis cf high 


speed phctccraphs to 


determine 


carbon icn velocities cf 


6 

2X1C cm/sec. & 


3C MH 


Q-switched 


laser was 


used for this 



exferiioent «ith a skewed pulse shape rising to its peak 
value in 10 rsec and ending at about 45 nsec. 

antcncv £2] established the ion velocity of a laser 

6 

picduced aluainum plasma to be between v = 5.2X10 cm/sec 

avg 

6 

and V = 8.6X10 cm/sec. In this wcrk a cylindrical 

max 

capacitor prcte placed 90 cm from the target was used for 

detecticc. The laser used was a Q-switched Nd:glass laser 
with a full width half maximum pulse length cf 25 nsec and 

1 2 

an energj cf 1.6J. The laser flux density was 1.1X10 
2 

S/cm . 
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Hei, Kelson 


, and 


Hall 


analysis cf heron 


nitride 


(BN) 


soca-lime glass 


targets 


arr 


sicnificartl j lowe 


r than 


oth 


foric; 







latle 2. Particle 



MOST PBC2ABLE 
VELOCITY 

s 

EM 2.25X10 



s 

A1 C 1.25X10 

a 3 

Scca-Line 

5 

Glass 0.73X10 



[59], through spectrcgtaphic 
, aluminum oxide (A1 C ) , and 

2 3 

ived at velocities that were 



er experimenters. 


In table 


Velocities in cm/sec 




AVEBAGE 


MJSXIflOM 


VELOCITY 


VEICCITY 


5 


s 


2.53X10 


2.76X10 


s 


s 


1.41X10 


1 . 53X10 


s 


s 


0.82X10 


0.69X10 



For this eaperiment a CO^ laser was 
nicrcseccrd pulse, 15 J. of energy and a re 

7 2 

density cf 10 S/cm . The microdensitom 

placed 4 cm from the laser-target impact pc 
for the very lew ion velocities in this 

the lew flux density of the CO^ laser used. 

Easev, et a-1. [7] measured the plasm 
carfccn iens *ith an electrical diagnostic p 
cylinder. Their results were plasma veloci 

7 

of 10 cm/sec. They also noted the velccit 

5 

■ region ct the plasma plume to be 5X10 c 



used 


with 


a 25 


latively Ic 


Vi flux 


et-er 


detect 


cr was 


int . 


The 


reason 


experiment 


tray be 


a ion 


velcc 


ity of 


refce 


in a F 


araday 


ties 


of the 


order 


y cf 


the 


opaque 


m/sec 


The 


laser 
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used was a Q-suitched Ndrglass laser with a pulse length 

12 2 

less than 15 nsec and flux density of 10 W/cm . Using 

ancther detection technique, Basov, et al. [6] again 
provided ion velccities for the laser induced plasna. The 
method of jlasma diagnostics used in this experiment was 
rather ingenious. Basov measured the optical density cf the 
plasma with a second laser at right angles to the first, 
several centimeteis out from the target surface. The 
teginning cf the plasma plume was detected by shining the 
second laser light through the plume. The opaque region of 

6 

the plasma plume was found to have a velocity of 6.3X10 
cm/sec fcr the aluminum plasma and the fast aluminum ions 

7 

had velocities cf the order of 1.6X10 cm/sec. A Q-switched 
Ndrglass laser was used with 150 Mli, 25 nsec pulse length. 

Dees any theoretical development also predict these 
velccities in the experimentally observed magnitudes? B. 
Steverdinc [!2] in the theoretical development of target 

V 2 

shock theory for incident laser flux greater than 10 K/cm 
shewed by tiecretical calculations that the tlcwoff velocity 

5 

is about HX1C cm/sec. He used a set of prcblem parameters 

similar tc these in the actual experimental situation used 
by Wei and Eall £5S]. An idealized laser pulse at 10.6 
microns kith an initial spike cf one microsecond duration 

7 2 

and maximum intensity of 2X10 W/cm was used in the 
calculations. The other parameters were target density 

3 

equal tc 5 gm/cm and the s-peed cf sourd in the target 

3 

material cf 2X>0 cm/sec. The problem solution was obtained 
using the cne'^dimensional theory of shock generation and 
plasma fermatien tc arrive at plasma velocities of the order 
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5 

c£ 4X10 cii/£€c. Ihe experimental results in the actual 

5 

ph 5 sical prctlem were 4.8X10 cm/sec. For air pressures 

telcw 0.1 ato. a blast wave will form at laser impact rather 
than a detcraticn wave. 
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III. theory 



In laser-target interacticn the question arises, wkat is 
the oiecharisD cr aechanisms cf plasaia production? Many 
facets cf the answer are as yet unexplained because cf the 
complexity cf the differential equations used to formulate a 
model cf tie problem. Laser radiation absorption ky the 
surface cf a solid metal target causes heat tc be deposited 
in a thir surface layer of the target material. The surface 
temperature ircreases and heat waves are propagated it the 
target. Ike keat cannot diffuse cut of the absorbing target 
layer as fast as it is being injected so the temperature 
continues tc increase until the target material begins to 
vaporize . 

The transition cf a gas into the plasma state involves 
various particle interaction processes. Collisions of the 
particles among themselves and interactions with radiation 
encompass most cf the energy transfer [20]. lonizaticr, the 
stripping cf electrons from atoms and molecules, is 
acccmplisled by the energetic (hot) electrons collidine with 
atoms and molecules. The inverse of this process is 
recombinatior . The atom or molecule has tc rid itself of 
energy. Three body recombination, in which two electrons 
collide simultaneo usly with an ion, one carrying away the 
excess energy and the other attaching itself, takes place in 
a cense plasma. 

There are three basic methods of radiation emission 
from particle interaction in a plasma; discrete radiation, 
recembinatien radiation, and bremsst rahlung . Discrete 
radiation arises from electron transitions from one energy 
level tc anctler in the same atom. Discrete radiation is of 
a single wave length, A complete picture cf this form of 
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xadiaticr cccsists cf numerous spectral lines. Another term 
used for discrete radiation is bound-bound radiation. 
Reccnibinatic r radiation is emitted when a free electron is 
captured by an ion. A lower energy state is assumed by the 
electron and the energy is emitted in the fcrm of a phcton. 
Eeccmbinaticn radiation is sometimes referred tc as 
free-bcutd iransiticns. Bremsstrahlung is radiation from 
interaction between a free electron and an icn in which the 
electron is only decelerated, not captured. Ihis is 
referred tc as free-free radiation. 

At the time of vaporization, the surface temperature 
of the target material will begin to depend cn the rate of 
vaporization; that is, the evaporation mechanism. Thermal 
diffusion continues tc increase the temperature of the dense 
sclid and ]iguid states still at the target surface but 
ceases tc play a significant role in the plasma temperature 
inc rease . 

In order tc develop the condition for plasma 
prcducticn, sufficient radiation energy density must be 
delivered tc the target. 

The electrostatic probe and associated theory was 
originally developed by Langmuir [31,36,55,56] and has 
evolved through a number of authors, Johnson and baiter 
[26], Eaksht £4], Lindberg [33], Virmont [56], Eyni [18], 
Chen [25], Swift [53], and Chung [14,15], tc its present use 
in plasma diagnostics, from space probes to fusion reactor 
research. The electrostatic probe has been used as a 
fundamental diagnostic tool for measuring plasma 
characteristics. The original work in its use and theory 
was dene by Langmuir (1924-1929) so the electrostatic 
probes, particularly the single filament probes are called 
Langmuir prehes. 

The electrostatic prehes themselves are relatively 
sinple devices, but the theory supporting the probe response 
is particulcrly complex (see Chen [25]). The electrostatic 
prebe basically consists of one or more small, metallic 
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electrodes (whiskers) that are inserted into the flasma to 
te characterized. The electrodes may be cylindrical or 
spherical ir shape, or may just be a plate or seme ether 
regular cr irregular geometric shape which suits the prebe 
suriace area calculation model and does net interfere with 
the plasoa flew or the plasma energy fields. The basic 
assumptions associated with electrostatic probe use, as 
listed b 5 Swift [53], are: 

(1) flectron and ion concentrations in the plasma are 

egual . 

(2) Electron and ion mean free paths are much greater 
than the electrostatic probe radius for a cylindrical or 
spherically shaped probe. 

(3) Electron temperature in the plasma is much greater 
than the icn temperature. 

(4) The probe radius is much greater than the Debye 
let gth. 

(5) The electron and icn velocity distributions are 
ha:twellian. 



In the single probe configuration, the probe operates 
by a single electrode inserted into the plasma and attached 
to a variable power supply . The power supply can be biased 
at various potentials, positive or negative with respect no 
the plasoa fci optimal probe signal response. Current at 
the probe is measured as a function of the prebe potential. 
The ether electrode (usually a plate) is a ground for the 
plasma. lie plate is normally fixed at a conducting section 
of the wall of the plasma confinement vessel. In electron 
discharge tubes the customary experimental arrangement uses 
the anede of the tube as the reference electrode. The 
Langmuir probe must have a continuous conducting medium from 
the plate tc the electrode which has distinct disadvantages 
when studjinc laser produced plasmas. 
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Ih€ ccuble probe configuration for plasma 
in vestigaticE is typically two symmetric electrodes that are 
in contact with the plasma. The double prcbe was first used 
by Johnscn and Walter [28], The current flowing in the 
plasma between the two electrodes is measured as a function 
of the vcltage applied between them. The double probe has 
the advantage ever the single probe in that there is nc net 
charge drair frem the plasma since the twe electrodes and 
the power supply of the double probe system form an isolated 
closed circuit, ie, , it is not grounded to any other 
equipment. Eoth electrodes are always negative with respect 
to the plasna and the maximum current drain by the prcbe on 
the plasma is limited by the ion saturation current tc the 
surface area cf the two probe whiskers. 

The theory cf electrostatic probes is complicated 
because the probe electrodes are themselves plasma 
boundaries. 2t and near the plasma boundaries the plasma 
characterist ics are changing, usually very markedlj. In 
other words, it is a grand dilemma. The necessity of 
getting inside cf the plasma with a prcbe that is by 
definition outside. Charge neutrality does not held near 
the plasna boundary. A thin layer called the Debye sheath 
[13] exists at the boundary where electron and ion rumber 
density differ from the values within the plasma and E 
fields are present. For a probe at plasma potential the 
sheath vacisles. 

Cher [25] and Swift [53] have summarized the 
theoretical results cf electrostatic probes in a systematic 
marnet- The experimental considerations when using 
electrostatic probes follow. 



|1) The surface of the probe whiskers should be clear 
of all centaninants. If the plasma temperature is lower 
than the wetk function of the probe whiskers, then the 
variation cf the work function over the prcbe surface will 
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affect the prcke characteristics. It is gccd procedure to 
ootgas the probe whiskers before starting an exper inental 
cata run, if enough current is available tc the circuit, by 
increasing the prcbe current enough to heat the whiskers tc 
incandescence . 

(2) Seccndary emission from the prcbe whiskers and 
arcing across the whiskers can also be a problem. The 
effects cf secondary emission are difficult to correct and 
can be avoided by using probe whiskers that have low 
emission ccefficients and by using lew probe circuit 
vcltages . 

(3> In a weakly ionized plasma, the presence of an 
electrostatic probe disturbs the nature of the plasma by 
lowering the plasma density in the vicinity of the probe 
whiskers outside of the Debye sheath around the prcbe. In 
high temperature pl-asmas, the probe can dissolve in the 
plasma, adding to the impurity problems in the plasma 
vessel. 

(4) Using a shielded probe for plasma diagnostics 
introduces still another boundary tc the plasma and all of 
the associated boundary problems. If using a shielded probe 
can be avoided while still satisfying the experiment design 
parameters, it should be. 

(5) The plasma probe signal is a function cf the 
surface area cf the probe exposed to the plasma. If the 
plasma energy is great enough to dissolve the probe, the 
surface area reduction of the prcbe must also be taken into 
account in the experimental design. 

(6) If the response characteristics of the prcbe are 
tcc slow tc follow the oscillations set up in the plasma, 
the probe signal will be useless. An approximation cf the 
expected plasma oscillation freguency should be made before 
designing the freguency response cf the probe circuit. 

(7) Election reflection, that is electrons that arrive 
at cne cf the probe whiskers which reflect back intc the 
plasma, can be a problem for accurate measurements in a 



29 



i' 

• - 
- \ - 

it i 

t • :■ .,: ■/ ^ »t , 

^ i >■ «V 

- • - a J '3 I t 

«*» ;.. • 

i - / ' • .i • . -. 

^ ■ ■: 

t % 



I 

..i It 

“I 




vr»- 



JJ I3>5J 

il 7 

. , aii^ o?s 

:BlftA»3>{C 

(SI 

?i>$ 9<it!>a; 
) ■» a3Dt:|j 

'iOBB »jt Ufr* 

o noi-aaisi 
ar 4£| 

't pai3»H9: 

^- U3 £3tM«ijlj 

.X»«s»i 

-.13 (HI ’j 

«»9u603tni 
^laoaa# »ifi 

.ai» 3 M 6 lA^ 

'«• j 

i}» 

aaiali 

4 If. 99Al3&t| 

( 

OI * j 

Ml 

• 3 »j 9. »iJ 

:« 39 t|n j 

7) ttfiaa! 



V 






«! iCi 



terucus flasia. 

(8) Ehctoelectron emissicn by the probe must be 
accounted fcr in experiments with a very tenuous plasma. 
When the plasma ion currents are small, photoelectron 
emissioE can have a dramatic effect on the picbe signal. 

(9) Ir small plasma discharges, the length cf the 
prcfce and the plasaa perturbations induced ty the prcbe may 
he large ccnpared with the length cf the gradient cf the 
plasma icn nucter density. The resulting effect cr the 
prche current must be compensated for. 

(lC-12) negative ions cause the collisicnless plasma 
theory assumptions to be invalidated, metastable atoms 
literate elections that alter the probe signal, and ions 
crtiting arcund the negative whisker cf the probe all cause 
prche signal ancmalies that must he realized in accurate 
plasma research. 

(11) The prote circuit, external to the plasma vessel, 
is a part cf the experimental considerations that is 
scnetimes cverloCked. Stray capacitance cf the 
electrostatic prche leads or in the circuit wiring is of 
great importance when attempting to detect high freguency 
signals. i floating double probe system will net follow 
rapid changes in the plasma potential unless the entire 
system has a small leakage capacitance tc ground. Ground 
Iccps and stray r-f signals are a particular problem arcund 
the prcbe circuit. 

The double prcbe theory used in the research 
discussicts can be found in Appendix £. A detailed 
fcrmulaticn cf the plasma temperature and density, the 
characteristic Debye length for the plasma, and the ion 
saturation current are all included in the appendix. 
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IV, EXPEBIMENIAL EESIGll 



In this research, the velocity cf a laser-prcduced 
Eiasma frcu an aluninuin disc target was investigated. It 
was a ccntinuaticn of research guesticns posed ty Ercoks 
[1C] and Schfcirzke and Cooper [^6,47], Ihe target was 



- 6 



suspended in 

-.6 

1,6X10 tcrr 



vacuum 



chamber at a pressure of 3X1C 

-6 



to 



being 


taken at 


1 


.8X10 




t 


err. 


sec (see Appen 


di 


X A) t 


c p 


ro 


vide 


laser 


pulse. 


Th 


e aver 


age 


1 
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that 


wer 
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ccities . 












e re: 
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€n 


the 


gnal. 


two-thir 


ds 


the 


spe 


€d 


of 


erist: 


ic speed 


cf 


the 6 


lec 


tronic 



in deteririnicg the plasma velocities. 

Ic calculate the time 
trigger, the pulse and the sigi 
light was used as the characteJ 
signal in the cables and the speed of light was used fcr the 
laser light ir the dehumidified atmosphere cf the laboratory 
air. The physical arrangement of the laser light path, the 

cillcsccpe are given in Figure 1. 
ns of the time relationship cf the 
ciggering circuit gave the following: 
path (AFGHJ) , the transit time cf a 
21.5 nsec. This was calculated using 

1 0 



delay ca 


ties 


and the 




Sa 


aple 


calcu lat 


plas 


ma 


sign 


al tc the 


Fcr 


the 


trig 


ger signa 


sign 


al 


at 


A was 2 


t = 


d/v. 


ie . 


11 



cm. Go = 350 cm, c = 3X10 

(the speed cf light in a vacuum). 

to to 

AG/c + GJ/(2/3)c = 120/3X10 + 350/2X10 

= 21.5 nsec 



cm/sec 
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FIGURE 1. THE TIME DELAY CABLES 
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Siailarlj, the other path 
i£ = 60 cii ano C£ = 300 cm. 



(AB+CDE) was calculated where 



10 10 

AE/c + CE/(2/3)c = 60/3X10 + 300/2X10 

= 17 nsec 



Note that tie distance BC does not appear in the calculation 
because that is the unknown velocity that was being 
measured . 

When the laser team struck the aluminum target at an 
ancle of 30°, a crater was formed. Figure 2 is a graphical 
representation cf the resulting hole. 



The 



Table 3. Crater Di 


mensicDS 






Dimecsicns 


10 Shots 


25 Shots 






A 


0. 160 cm 


0.175 cm 






E 


O.hIC cm 


0.4C6 cm 






C 


0.105 cm 


0.115 cm 






c 


0 . 1 GO cm 


0.225 cm 






E 


0.704 cm 


0.701 cm 






meascrements in 


Table 3 were 


made with a 


caliper 


and 


rcmeter. The laser 


pulse removed 


0.01 cm cf 


the hcttom 


the crater hole on 


each shot (di 


mension E divided hy 


the 


ter of shots) . Tc 


calculate the 


laser flux 


density 


at 



the target, it was assumed that the pulse length in seconds 
(t) egualed the half width of the pulse (25 nsec) and the 
diameter (d) cf the beam incident to the target was 
dimension C from figure 2. laser energy equaled 8.6 J on 
the average, sc E = 3.6 J, t = 25 nsec, d = C.110 cm. 

2 

(0.055) 



2 -9 

fluj density = E/t tt (d/2) = 8.6/25X1C 



10 2 

= 3.62X10 W/cm 
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FIGURE 2. 



THE LASER-TARGET 



INTERACTION CRATER 
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The possibility of drilling completely through the alcirinuni 
disc had tc be guarded against so that all of the laser 
prcduced plasma came out the front of the target. The 
thickness of the disc measured with a micr cmeter was 0.635 
cm (0.25” machined aluminum). If s was the depth the 
crater must reach to drill a hole all the way through the 
target and each laser shot removed approximately 0.01 cm of 
the target material, then it was estimated that 70 to 75 
shots would make a hole through the disc (see Figure 3). 

s = d/cos 30° 

= 0.635/0.866 
= 0.733 cm 
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vacuum chamber (see 
• was inserted in a probe 
:ectangular brass plate. 

1 top of the chamber and 
a bench mark about in 
prcbe. Ihe plate was then moved so the prc 

:e. A paper cent! 
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brass 
and the 



plate could 
probe dist 



from the surface 
arrangement was 
Figure 4) . The 
holder which was 
An aluminum butt 
the brass plate 
the center cf the 
be whiskers just 
meter graph scale 
ed down cu the 
then be moved a 
ance would be 
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LASER TARGET IMPACT ANGLE 
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COAXIAL CABLE TO 
THE DC POWER SUPPLY 




FIGURE 4. THE PROBE DISTANCE SCALE 
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Ihe pcver supply used for the electrostatic probe 
circuit (see Appendix B for prcbe construction, theory and 
circuit desicE) was a Hewlett Packard 6218A power supply 
with a naxiiiuai voltage of 50 volts and maxinium current of 
2C0 mA. iarly in the research using the electrostatic 
prcbe, it was determined that it was desirable tc ccmpare 
the signal peak values for different values cf probe bias. 
A relative nuiber equal to the measured peak voltage divided 
by the bias vcltage could then be used for data reduction 
and graphing. The prcbe response, as graphed by Brccks flO] 
in previcus researc-h using the same probe, indicated the 
recicn frcrc -10 volts to -20 volts DC bias to have the 
linear characteristics desired for accurate compariscns (see 
Figure 5 . ) . 

The pcints on this graph were an average of several 
hurdred laser data shots, especially at -IS volts acd -20 
vclts. Ihe standard deviation cf the main plasma peak value 
fcr those bias vcltages was abcut 4 volts in each case. The 
average values did match up to Brooks* [10] graph of the 
sane characteristics hut no mention was made by Brocks of 
the accuracy cf the data used tc plot the curve in Figure 5. 

When the probe bias vcltage was determined, the data 
pictures could be taken. Five oscilloscope time scales (x 
axis) were used tc take the data from the probe circuit, 20 
nsec/cm, 50 csec/cii, 100 nsec/cm, 200 nsec/cm, and 500 
Esec/cm. Fcur vcltage scales (y axis) were used tc record 
the probe signal, 500 millivolts/cm, 1 vclt/cm, 2 vclts/cm, 
and 5 vclts/cm. Ihe 20 nsec/cra and 50 nsec/cm scales 
displayed the early plasma signals; that is, the fast acving 
particles that arrived prior tc the main plasma pulse. The 
IOC and 2C0 csec/ca scales showed the main plasma pulse and 
the 5CC nsec/cm scale shoved the late plasma signals, i.e. 
the particles with slew expansion velocities. 
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FIGURE 5 




THE PROBE RESPONSE WITH THE PROBE 2 CM FROM THE 
TARGET SURFACE 
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FIGURE 6. DATA PHOTO USING A 62 NSEC TIME DELAY CABLE 
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To te able to consistently record the signals in the 

saae time reference, a zero time mark, t , had tc be 

0 

displayed cr each photograph. The inverted laser pulse 

signal was used for the time reference. Initially, a long 
tine delay cable (62 nsec and 48 nsec cables were used) was 
added to the probe circuit so that the probe signal followed 
the laser pulse signal by the time length of the delay 
cable. This technique worked very well for the early 
signals (20 nsec/cm time scale) but was fcucd to add little 
tc the larger oscilloscope time scale selections. Cn a 50 
nsec/cm scale the delay cable cnly delayed the signal cne cm 
and when a SCO nsec/cm scale was selected, the 0.1 cm delay 
cn the cscillcsccpe graph was not even discerrable. 
Eventually, the following techniques for timing were 
established as the method of extracting data from each time 
scale. 

(1) The 20 nsec/cm time scale was easily read to 
4 nsec. The laser pulse covers about 50 nsec from start to 
finish, sc 48 nsec and 62 nsec delay lines were used and 
zerc time fcr the plasma signal was assumed to be at 5C nsec 
or 60 nsec respectively. Measurements were made starting 
2.5 cr 3 cm frcm the edge of the graph (Figure 6.). This 
scale can read all probe signal peaks to about 80 nsec. 

(2) For the 50 nsec/cm time scale, the best 

accuracy cf measurement was 10 nsec. Since the 20 nsec/cm 

scale could read the signals tc about 80 nsec, the 50 

nsec/cm time scale was used for signals frcm 80 to 500 nsec. 

Fcr this scale and all higher time scales, the time delay 

cable was removed and a replacement cable was used tc match 

the trigger pulse up with the probe signal (see Figure 1) . 

The peak cf the pulse signals was used as the reference time 

t where, 

0 
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40 + 5 = 45 nsec. 



t = 
0 



All of tb€ time values were aeasured from t , addina 45 nsec 

0 

to get the actual time of the peak signal. For example, 
suppose a prcbe signal peak was 115 nsec from the laser 
pulse peak. 



t =t+t =115+45 

peak 0 

= 150 nsec. 



(3) The ICO nsec/cm scale cculd he read with 20 
nsec accttacj. The data from the 100 nsec/cm time scale was 
read like the 50 nsec/cm scale. The 100 nsec/cm time scale 
was used for peaks located from about 150 to 800 nsec. 

(4) The 2C0 nsec/cm scale cculd re read with 40 

nsec accuracy. Now even the t value determined in (2) 

0 

above is heccning insignificant. However, because of the 

ease cf locating the peak of the laser pulse, the same 
method was used to determine time relationships (see figure 
7) . The 2C0 nsec/cm scale could be used for signals from 
20C nsec tc about 1600 nsec. 

(5) The 500 nsec/cm time scale was used tc see if 

there were any very late signals. It cculd detect signals 

out to roughly 2,500 nsec. The time values accuracy for the 

500 nsec/cm scale was 100 nsec. The t value was not 

0 

significant to these figures and the time value of the peaks 
was read directly from the laser pulse peak (see Figure 8) . 
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LASER 

PULSE 




(a) DISTANCE FROM THE TARGET SURFACE TO 
THE PROBE EQUAL TO TWO CM 




(b) 



FIGURE 7. 



DISTANCE FROM THE TARGET SURFACE 
TO THE PROBE EQUAL TO FIVE CM 



TIME REFERENCE ON A 200 NSEC/CM TIME SCALE 
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FIGURE 8. THE PLASMA PROBE SIGNAL ON A 500 NSEC/CM SCALE 
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To take data from the graphs, a magnifying glass and 
Plexiglas certimeter grid scale were used. The magnifying 
glass facilitated more accurate readings cf the clastic 
scale which cculd be lined up with the time reference pcint 
mere readily. It also helped to reduce the eye strain 
pursuant tc reading the graphical data. 

The plasma probe signal occasionally had a high 
freguency ncis§ signal superimposed on the laser pulse 
signal (see figure 9) . The noise signal only appeared 
during tte time cf the laser pulse. To determine what the 
signal was, the frequency was approximated as follows. From 
a 20 nsec/cm data shot, tvio oscillations were measured tc be 
C-ti cm long. That indicates the time constant, t , to be 

approximately 4 nsec. The freguency then is. 



- 9 

f = 1/t = V'^XIO 

s 

8 

= 2.5X10 cycles/sec. 

Unwanted electrical grounds were first suspected. Herrison 
[35] was consulted for grounding problems, but all were 
discounted. It was first thought that the noise was the 
electromagnetic field radiation from the Pcckels cell being 
picked up by the MgO photodiode detector (K-Dl) . The 
detector fed the laser pulse signal tc the oscilloscope. 
The noise seemed to be cn the laser pulse signal. 
Intrcducirg a Icng time delay line (62 nsec) or the 
electrostatic prebe signal separated the noise from the 
laser pulse sc the noise had to be coming from the probe. 
It was then conjectured that the probe vas acting as an 
antenna receiving a signal frem the Pockels cell, but that 
tcc was discounted because the signal disappeared when the 
prebe was pulled up away from the target. No signal was 
received when the laser was fired with the probe in that 



45 



i 



* ! t* 

■n 

I * 10 $ 

<i : 

•n^Up»7l 

' oa<t2n& 

V iBOI^il 

^en 01 ft 

• 0 !•) 




•>iaft«gq| 
stt ti€ 
.lusaB 
j a>tf» 

Jt#f 

»-(fr 

ia»( 

I* 4 | 



I 



pcsiticc, which contradicted the hypothesis that the probe 
was an antenna. Further investigation, to see how the noise 
signal changed as the laser burned a deeper hole intc the 
target, indicated that the noise signals had something tc do 
with the laser target interaction. For the first fifteen 
shots the noise signal was consistently present. Is the 
laser hcie get deeper, the signal disappeared. After the 
15th shot, the signal rarely appeared and after the 40th 
shet on a single target spet, there was never a noise 
signal. Alsc, from a highly polished aluminum target 
surface (almost mirrer quality), the noise signal was even 
mere marked than from a dull surface. The conclusicn was 
that the rcise signal was due to either X-ray emission from 
the laser-target interaction or some sort of diffuse 
reflection (twinkle) from the target surface. To minimize 
the noise, the target surface was ••cleaned” with an 
unfocused laser shot before each data run. The cleaning 
reduced the persistence of the high frequency noise signal 
but it did not totally remove it. 

The data collected from each photograph were tabled 
for data reduction calculations (see Appendix C) . For each 
picture, a unique 10 digit number was assigned tc identify 
it. The shet number per target spot, N , on which the 

C 

picture was taken was recorded. The distance of the probe 

frcoi the target surface and the probe bias voltage were 
recorded. Cc seme pictures, as many as three peaks could be 
dis tingui she c . The time and voltage of each signal peak was 
recorded. The laser energy was recorded and the error of 
measuremert cf the peak times and peak voltages were 
recorded . 
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FIGURE 9 




HIGH 

FREQUENCY 

NOISE 



THE HIGH FREQUENCY NOISE SIGNAL SUPERIMPOSED ON 
THE LASER PULSE SIGNAL 
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A linear time tc distance relaticnshif was assumed and 
calculations were made on the basis cf this assumfticn, 
ie. d = vt. The velocities agreed with previous work at the 
approximate distance of the probe from the target surface (1 
to 5 CB) . BcKee [44] determined the plasma velocity tc be 

7 

10 ciE/sec and Brooks [10] found laser produced plasma 

7 8 

velocities cf 5.SX10 and 1.1X10 cm/sec. Demtrcder and 

7 

Oantz [11] measured velocities of 10 cm/sec at 170 cm and 

6 

Hei, Nelson, and Hall [59] measured velocities of 10 cm/sec 

at 4 mm sc the linear velocity assumption was a valid one 

fcr the region 1 cm to- 5 cm from the aluminum target 
surface. 

The eguaticn used to calculate the velocity, v, from 
the probe signal data was v =■ d/t, where d is the distance 

the plasia had to travel and t is the tine of the plasma 

prcbe peak. A velocity was calculated for each probe signal 
peak. Ihe distance, d, was assumed to egual the distance 
from the surface of the aluminum target to the prcbe, d^, 

plus the distance from the bottom cf the hole tc the surface 

of the tarcet, d , where d = 0.01 N [cm] because each shot 

2 2s 

removed C.C1 cm of the target surface in the crater. If 1 

2 

were not taken into account, the deep holes (N > 40) would 

s 

have large measurement errors before the data reduction 

calculaticns are even begun. 

Error analysis of the photographic data calculeticns 
of velocity and distance was also necessary tc insure that 
the measured guantities were dependable, A simple scheme 
fcr error analysis found in reference [3] was used. The 
measurement errors for distance d were <5 d = ± 0.05 cm. 
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tiaie errors, dependent on which oscilloscope tiise scale 
used, were: 

lable 4. Cscilloscope Time Measurement Errors 



SCALE 

/ iO nsec/cm 

50 nsec/cm 
ICO nsec/cm 
200 nsec/cm 
500 nsec/cm 

velocity calculaticn error 
V = d/t, sc 



TIME EEEOR (dt) 

± 2 nsec 
± 5 nsec 
±10 nsec 
±20 nsec 
±50 nsec 

as . The velocity formula 



6 V 



l3v/9dl 6d + |3v/9tl 6t 

1 1/t 1 6d + 1 -d/t^ 1 6t 



ample oaloulaticn of the velocity error for velooities of 

6 

order of 2X10 cm/sec, with t = 50C nsec, d = 1 cm, 
= 50 nsec, and id = 0.05 cm gives: 



-9 2 - 9 

iv = (1/500) (0.C5) + (1/(500X10 ) ) (50X10 ) 

s s 

= 10 + 2X10 

5 

= 3X10 cm/sec. 
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Tatle 5. Results Velccity Error Calculations 



V [cm/sec] 


t 


d 


[cm] 


^t 


dd [cm] 


dv [cm/sec] 


<5v/v 


6 

2X10 


5C0 


nsec 


1 


50 nsec 


0.05 


5 

3.0X10 


155J 


6 

6X10 


33C 


nsec 


2 


10 nsec 


0.05 


s 

3.3X10 


5.555 


7 

5X10 


IOC 


nsec 


5 


5 nsec 


0.05 


6 

3.0X10 


655 



So the measurenent errors of the velocities were not very 
important ficm a statistical point or view. for data 
analysis, eticrs of up to 305^ were accepted as good data. 
Ihe total cumber of velocity data points was reduced from 
585, the numter of peaks on the data photos, to 171, the 
number that ultimately survived the error analysis gate. 
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V. 



RESCLTS 



Ihe exf eriioeDtal results of the research completed for 
this thesis are presented in this section. Each aspect of 
the results has been subdivided into a subsection fcr easy 
identification. The results of this thesis were extracted 
frca graphs chtained from computer analysis of the data. 
Pcur plasma ion velocities were found corsistently ir each 

6 6 

Shot series; 4.0X1C cm/sec, 5.6X1C cm/sec (the bulk cf the 

7 7 

plasma), 1.5X10 cm/sec, and 3.7X10 cm/sec, and a fifth, 
e 

10 cm/sec, appeared in a limited number of the data runs. 

A. C5ATIE EIKENSICNS 

Ihe crater dimensions were taken free two craters for 
measurements tc approximate the amount of target material 
removed hy each laser shot. Ihe results were tabled in 
Table 3 and a graphical representation of the crater is 
shewn in Figure 2. The deep crater had a cylindrical 
diameter cf C.11 cm. Each shot removed approximately 0.10 
cm cf target material. Ihe surface hole cf the crater 
differed ir dimensions slightly frem the deep, cylindrical 
shaped hole. It had a diameter of 0.17 cm. Ihe difference 
was probably due tc the surface boundary at the top cf the 
crater only being bonded to the target cn cne side while the 
deep crater material had target material cn three sides 
which reduced the erosion at the edges of the target spet. 
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B. PROBE SIGKJL CHANGE ASSOCIATED WITH THE EIFTEENIH lASEH 
SfiCI j 

There seemed to be a consistent probe signal 
difference associated with about the fifteenth laser shot 
number that changed the laser plasma characteristics when 
the laser-target interaction angle was 30°. It was detected 
in this research with an electrostatic probe, Shewchuk [49] 
detected it with X-ray photodiodes, and Hilliamscr [61] 
detected it with magnetic probes. 

To explain this anomally it was speculated that the 
30® laser team incidence angle was hiding the incidence 
impact area from the probes after the fifteenth shct. A 
guick calculation using the laser cratering statistics of 
Table 3 showed that by approximately the fifteenth shct the 
impact area is hidden. Additional support tc this 
hypothesis was lent by Shewchuk’s research. With the X-ray 
prches in a position normal to the target surface, the X-ray 
signal began to get weak after the fifteenth shct. One 
X-ray prche was continually collecting a greater signal 
strength than the other. When the probes were pcsiticned 
nearly parallel to the laser beam so that they could nearly 
lock directly down the crater there was no signal charge no 
matter hov many shots were taken on a single crater. The 
consistent X-ray signal indicated that there was nc special 
significance tc the fifteenth laser pulse ether than the 
mechanical bJccking of the signal by the crater edge. 
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c. 



PLASHi ElECIECK lEMPEEATUEE 



ElectrcD plasma temperature can be found fron the 
slope of the electrostatic prcte characteristics (see Figure 
5) at the graphical origin, (0,0) , using an eguatior. for 
electrostatic probe diagnostics from Chen [25]. 

dl/dTI = (e/kT ) (i i )/(i + i ) 

e »+ 2 + 1 + 2 + 

= ei/2kT 

e 

where i = i = i because the double probe used was 

2 + 

syumetric. Ihe measured value of dl/dV; 

dl/dV = 1.5/5 = 0.3 amps/vclt 

by extrapolation of the graph in Figure 5 and comparison 

with Brocks' graph [10] of the same data. Ihe value of i is 

i = BV where V , the probe saturation voltage, is equal to 
s s ■ 

30 vclts and E, the signal resistance, is equal to 1 ohm. 



T 

e 



-1 

(ei/2k) (dl/dV) 



- 1 

= (11605/2)30(0.3) 

4 

= 58X10 OR 
= 48 eV. 



Sc, the plasna electron temperature is approximately 50 eV 
two cm frcm the target surface. 
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D. PLASHl EENSIIY RESULTS 



The plasffa 
plasma velccity 



- 3 

density, n [cm ], for each characteristic 
i 

is derived in Appendix B frcm the equation: 



n = I /Ze (1 + s) A V . 
i + p LP 



Appendix B defi nes all of the terms c£ the equation 
and makes the necessary approximations to arrive at the 
density fcr each prcbe signal velocity (see Figure 13) . The 
analytic ccrrespcndences of velocity to density were found 



< 13 -3 6 
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CCX1C cm/sec to 


4.09X10 


cm 


, 5.5 


7X10 cm/sec to 


*/ .34X10 


1 3 


-3 7 

cm , 1.52X10 


cm/sec 


to 


1 3 

10 


-3 7 

cm , 3.11X10 


cm/sec 


tc 


12 -3 

1.C3X10 cm , 


8 

and 1.C1X10 


cm/sec 


1 1 

to 1.C6X10 



-3 

cm 
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E. 



VELOCITY DISTRIEOIIOIJ ANALYSIS 0? THE EBOBE SIGNALS 



The velocity distributicn of the plasia particles is 
an insportant consideration in containnent scheces for 
plasmas. The container wail interaction with the plasma can 
te predicted if a velocity distribution model for the plasma 
is accurate. The point of this section is to analyse the 
probe temporal signals for velocity distribution. 

In order to deduce velocity distr ihutions frcn the 
prche current, • the time function must be analyzed. If the 

particle velocity is given by v = d/t, then v is 

z 

prcpcrticral to 1/t . The one ditnensicnal Maxwellian 



velocity distribution is: 

1/2 2 2 

dN/dv = p (V ) = N(l/v { n} )exp(-v /v ) 

X V X m X m 



V2 

where v = (2hT/m) 

m 

N = number of particles 
k = Boltzmann's constant 
T = plasma temperature 
m = particle mass. 

A Gaussian time distribution, f(t), is of the form: 



where 



2 2 

f(t) = A exp(-’(t-t ) /2s ) 

0 0 

A = the peak value 
0 

t = value of t at the peak 
0 



and 



s 



standard deviation. 
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FIGURE 10. 



DATA PHOTO FOR GRAPH OF 



ION VELOCITY DISTRIBUTION 
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The staccard deviation is a parameter vhich contro 
relative flatness cf the bell shaped curve. A small 
for s results in a sharp peak and a large value of s r 
in a broad, flat shape. 

Nov if all particles of a velocity grcu 

considered tc start from a single point at cne time, 
the Maxvell Eclt2mann distribution of velocities 
direction, then the distribution of the particles 
respect tc arrival times at the probe can he deduced a 

dt/dt = (dN/dv ) (dv /dt) , 

X X 

and since v is related to the arrival time thrcu 

X 

relationship v = d/t, then 

X 



dv /dt = -d/t , 

X 

2 - 1/2 22 2 

and dS/dt = {N ( jt v ) exp[-(d/t) /v ]) (-d/t ) 

m m 

-2 2 2 
= (constant) t exp[-(d/t) /v ]. 

. m 

This fcrmula does not analytically approximat 

irent limit. In the course o 
>e signals, however, appear 
[see Figure 10) . Further gra 
)s indicated that velocities 
reasonably fit Gaussian curve 
:he- photo in Figure 10 was g 
curve (A-B) successf ully . C 
the photo. In Figure 1 
out 





Ih 


is fcrmula 


Gaus 


Sian 


f era 


in 


rese 


arch 


the pi 


cbe 


appr 


each 


a Gaus 


sia 


analysis 


cf the 


da 


f ten 


th 


e p ich 




Figure 1 


1, the 


dat 


and 


fit 


ted tc 


a G 


was 


the 


actual 


cur 



results cf subtracting 



the particles with velc 



that belcncec to the next peak (curve E) made the Ga 



fit ever tetter, 
eguation was; 



For the curve A-B the Gaussi 



Is the 
value 
esults 

p are 
with 
in one 
with 



ch the 



€ the 
f this 
ed to 
pbical 
taken 
E . In 
raphed 
cr ve A 
1 the 
cities 
cssian 
an fit 



57 



•/I Jilai 

..-2 z 

ad £ tft 
:K 

til si# 
^ issziJi 
#:>»|ss.x 



'll* &s« 

ttoi#sE*T 



I 



ifr 

ias&iss; 

u3S31|jg 

tf «»$l 

-*ti 

tJlSt 

«i i«ia 

S- ffs 



( 



z 

f(t) = 2. 53exp (- {t-150} /4419) 
arc for curve E: 

2 

f(t) = 3. lOexp (- {t-160} /3804) . 

Fui;tb€r acalisis tc explain these unexpected results is 
required . 
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F, PLASMi PiSBIICLE VELOCITY GSAPHICAL RESDIIS 

Ih€ tctal picture of all cf the plasma velocity data 
with measurenent errors less than 30% was graphed in Figure 
12. Each cata peak was used tc calculate a velocity and 
each signal peak value, V, was divided by the DC prcte bias 
value, V , tc state each signal value as a unitless 

rc 

fraction, V/V , for inclusion of a range cf values for V 

IC DC 

There is a marked visual effect of the distribution of 

points it Figure 12 that immediately suggests that there are 

several discrete characteristic plasma particle velocities. 

Figure 13 carries the analysis one step further, shewing the 

graphical results of averaging all of the velocities and 

peak values in groups around each particularly dense region 



12. 


6 






slow signals (less than 


5X10 ci/sec) 


did 


not 


s valid points when 


the naximum alio 


wabl e 


data 


30j!. If the error was allowed to be 50% 


then 


the 


als remained in the 


analysis. That 


is 


what 


be slow signal velocity 


in Figure 13. 


Particle 



velocities were found to be 4,0X10 ± 1.2X10 cm/sec with a 

6 6 

signal peak value (V/7 ) of 0.21 ± 0.10; 5.6X10 ± C. 25X10 

DC 

cm/sec with a signal peak value (V/V ) egual to C.50 ± 

DC 

7 7 

C.13; 1.5X1C ± 0.26X10 cm/sec with a signal peak value 

7 7 

equal tc C.19 ± 0.08; 3.7X10 ± 0.76X1C cm/sec with a 



8 6 

signal peak value egual to 0-05 ± 0.04 and 1.0X10 ± 0.5X10 
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The 



cm/ssc with a sicnal peak value equal to 0.03 ± 0.03 

a 

very fast (1.CX10 cm/sec) plasma velocity was taken f 

data sample cf cnly three points out of 171 data p 
Even tbouch the error bars for that data in Figure 1 
small, the prctable error is not, due to the paucity o 
points. 



In Appendix B the 



ion density, n , was table 
i 



function cf the 



plasma particle velocities. 



V . 

IP 



Ecltzmann egcation 
1(f) = 

If E = 

then n = 



2 

where v = 

m 



is : 



-E/kl 

n e 
0 



IDV /2 

2 

n^exp (- {m V /2kT} ) 
2 2 

n exp(-v /V ) 

0 n 



m/2kT. 



6 1 



torn a 

oints. 
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f data 

0 as a 
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FIGURE 14, GRAPH OF ln(n.) VERSUS 
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FIGURE 15. GRAPH OF ln(ln(n)) VERSUS ln(v 
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a 2 

So, lt(n) = In (n )-(v /v ) 

0 01 

or a linear eguaticn of the form y = ax+b 

where y = In (n) 

b = In (n ) 

0 

2 

a = -1/v 

m 

2 

and X = V . 

2 

If the graph cf In (n) versus v yields a straight 

I 

lice ther the set cf points in Figures 12 and 13 fit a 
Ecltzmanc distribution. 



2 

latle 6. Data Points for ln(n) vs. v 



2 



V [cir/sec] 


V 


In (n) 




6 

5.57X1C 


13 

3.1X10 


31.93 . 




7 

1.52X1C 


1 4 

2.3X1C 


29.93 




7 

3.7121C 


IS 

1.4X10 


27.66 




3 

1.C 12 10 


1 6 

10 


25.41 




The graph 


2 

of In (n) vs V 


, shown in figure 14 did 


not 


produce a straight 


line, so the 


2 

relationship of n to v 

i 


is 


not a Ecltzaann 


distribution. 


A plot cf In (In ( r) ) 


vs 



2 

ln(v ), however, did produce a reasonably straight line (see 

Figure 15). Sc, the theoretical relationship is an 
exponential power relation. 
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If 




1 


c (n) 


= C (V ) , 


where 


c and 


K are co 


nstants 


m 










2 


K 




2 






Sc, 


In (In 


(n)) 


= In (c (V 


) ) = 


In (c) + 


Kin (V ) , 


is 




the 


graph 




of a 


straight 


line , 


y = ax 


+b, with 


y = In ( 


In (n) ) 


and X 


= It 


(V 


2 

) . 


The graph 


would 


then 


give a 


gecme 


trical 


value 


for 




F egual to 


the si 


ope (s 


ee Figure 


15) . 


If- the 


slope 


is K 


n 


ct egual to 1, 


then 


the B 


oltzmann 


eguati 


on is 


not a 


sclu 


ti 


cr . 















Ih€ crafh of ln(ln(n)) versus ln(v ) is a straight 

-2 

lice with slope egual to -3-86X1C . So, the relationship 

between n and v is an exponential power relationship, rot a 
Ecltzmann distribution as postulated. 

Error analysis- indicates that the probable error in 
the y ooorcinate of Figure 15 is negligible and the error 
tats in tie x direction (about 3551) are such, that if the 
data point is taken to be at either extreme, it is still on 
the line in figure 15, indicating a very reliable power 

2 

relationship between ln(n) and v . Graphical analysis 
yields : 

2 - 0.0386 

In (n) = c (V ) 

2 - 0.0386 

and the number density is proportional to exp(c(v } 

On the assumption that all ions detected are aluminum, 
the kinetic energy differences between the characteristic 
velocity groups can be compared with the ionization energies 
of the multiply ionized states of aluminum (shown in Table 

7 6 

8). For example with v = 1.52X10 cm/sec and v = 5.6X10 

2 i 

ca/sec, two of the characteristic plasma particle 
velocities, the resultant energy difference for interval I 
on Figure 14 is: 
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wh€I€ IS i£ 



acc 



£C 



z z 

E -I = (mv /2) - (inv /2) 

< X 2 1 

tt€ mass of an aluminum atom, 

-24 

m = 27(1 . 67X 10 ) gm 

2 2 14 2 2 

V -V = 2X10 cm /sec 
2 1 

-9 

£ -E = 4.53X10 erg 
2 1 

= 2.71 KeV. 



Similarly E -E and E -E were calculated for intervals II 

3 2 4 3 



and 


III. Ecr an 


aluminum 


mass E “E = 

3 2 


16 


KeV 


and 


E -E 

4 3 


= 122.7 iCeV 


If m is 


assumed to be a 


hydro 


gen 


mass. 


then 


E -E = 590 eV 

3 2 


and E -E 

4 3 


= 4.54 KeV* If 


the 


part 


ides 


were 


assumed to be 


electrons E 


-E = 2.51 eV. 









4 3 

The value of the change in kinetic energy resolved 
from these calculations is similar to the difference between 

12+,13+ 11+ 

the icnizaticE energies of Al and Al ions. The 



icnizaticr energies are listed here for comparison (all of 
the ionization energies are rounded to the nearest whole 
number) . 



Table 8. 



Aluminum Ionization 



Energies 



♦ 



Al 


6 


eV 


2 + 

Al 


19 


eV 


3 + 

Al 


28 


eV 


4 + 

Al 


120 


eV 
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154 eV 



s + 

A1 

fil 190 eV 

7 + 

A1 240 eV 

8 + 

A1 285 eV 

» + 

A1 330 eV 

10 + 

A1 400 eV 

1 1 + 

A1 440 eV 

»2 + 

A1 2090 eV 

13 + 

A1 2300 eV 



It is DCt clear why the velocities are grouped, tut 
locking at the aluninum ionization energies, four groups can 

+ 2 +^ 3 + 

te readilj approximated; Al in group one, Al in group 

4+-11+ 12+^13+ 

two, Al in group three, and Al in group four. 

If a local Boltzmann equilibrium is assumed to hold 
for each irterval in Figure 14, then the kinetic energy of 
2.7 KeV derived above can be equated to a temperature. Ihis 
temperature does not correspond with the X-ray energy of 
several hundred electron volts detected by Shewchuk[ 49 * , but 

+ 13 

it does correspond very well to the Al ionization energy. 
The possitility exists that the velocity groups result from 

I 2+,13 + 

the degree of electron stripping. For example, Al 

8 

ions correspond to a signal velocity of 10 cm/sec, group 

7 

three icrs correspond to velocities of 3.7X10 cm/sec, group 

7 + 

tkc ions correspond to velocities of 1.5X10 cm/sec and Al 

6 

ions ccrrespotd to 5.6X10 cm/sec. 
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Th€ target material should be changed from alunirum to 
seme other netal element whose electron separation energies 
are also known. This would either lend further support to 
the group energy hypothesis or cause the hypothesis to be 
rejected. Pirther study of this phenomenon is required for 
a complete understanding of the theory. 



G. CIHEB VEICCIIY BESDLTS 



There seems to be a mathematical pattern associated 
with the four plasma velocities. For convenience of 
nctation let the characteristic plasma velocities be 
represented ky: 

6 

a = 5.6X10 cm/sec 

7 

b = 1.5X10 cm/sec 



7 

c = 3.1X10 cm/sec 



and 



8 

d = 10 cm/sec. 



Then, fractionally comparing each velocity tc the 
velocity closest tc it in magnitude with the faster velocity 
in the numerator and the slower velocity in the denominator, 
the relationships ate; 

c/c = 2.7 
c/b =2-5 
h/a = 2.7. 

This relaticnship is a noteworthy constant. Its 
applicability, however, is not apparent. Similarly, other 
velocity relationships are also geometrically ordered; 

d/b = 6-7 
c/a - 6.6 
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There is considerable implication that these 

relaticns are significant. Their application to 
theoretical development is net clear. A pc 

explainaticn may lie in the aluminum ionization 
groups discussed previously. Pursuit of this conjectu 
not carried any farther in this thesis- 



number 
any 
ssible 
energy 
re was 



H. 



CCKPCSIIE 



EECEE BESEONSE 



The main plasma signal velocity and shape cempared 
well with previous work. However# the early signals which 
were thcrcuchly investigated were neither consistert nor 
reproducible. Iseccr [27] found that the early signals were 
strcngly dependent on the surface conditions and the 
ccmpcsiticn cf the gaseous atmosphere inside of the sacuum 

- 5 

chamber. Isenor was working in pressures from 10 terr to 

- 3 _ e - 9 

10 torr. Eerteus [41] found that at 10 to 10 torr 



there were nc early signals. Porteus indicated that the 
impurities it the space inside of the vacuum chamber could 
be the cause of the early signals; vaporized oil frem the 
diffusion pump# residual atmospheric air gases at the low 
concentr atiers inside of the vacuum chamber# aluminum vapor 
frcB previevs laser shots# aluminum oxide vapor or the 
target surface# and even grease vapor from human handling of 
the target while putting it in the vacuum chamber and taking 
it cut. All of these sources of contaminants tend to 
distort the probe signal from the aluminum plasma. 
Hydrogen# from a molecular layer of diffusion pump cil on 
the target surface, could also be the cause of the early 
signal. 

Frcff locking at many electrostatic prebe pictures, an 
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average electics tatic probe response signal seemed to 
suggest itself (see Figure 16). The peaks get larger when 
the probe was closer to the target surface and grew snaller 
as the prebe was moved farther away, indicating that the 
plasma spread out as it moved away frem the target surface. 
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FIGURE 16. AVERAGE ELECTROSTATIC PROBE RESPONSE SIGNAL 
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VI. CCNCLCSIONS 



In ccEclusiqn the results of the research represented 
bj this thesis are to be summarized. 

Ihe laser-ptcduced crater was made at the rate of C.01 
cm per laser shot. The crater cylinder diaireter was 0.11 
cm. At the fifteenth laser shot on a target crater, the 
impact surface was shielded ficm the probe when the laser 
incidence ancle was 30°. 

The plasma electron temperature was calculated 
analytically at two cm from the target surface to he about 
50 eV. 

11 13 

Ihe plasma icn density varied from 10 to 7.34X10 

- 3 

cm for the streaming plasma. These densities were 

calculated ficm data using a double electrostatic probe. 

Icc groups cf different velocities were observed. 
Within the main group, the plasma particle velocity 
distributicns were shown to be approximated by a Gaussian 
distribution. The icn number density in each velocity group 
was demoEstrcted to obey the Boltzmann equation for limited 
velocity ranees. 



z z 

r (v) = n exp[ -v /v ] 

0 m 
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Belatively large teoiperatures were calculated using 
Ecltznanr a^f roxinatiocs. An almninuoi ionization energy 
theory was stated to offer an explanation, tut more 
questions need to ie answered before a conclusive argument 
can te offered. 
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VII. FOBTHER BESEECH 



Ihe areas warranting further study that arose during the 
course ot expetinentation are: 

(1) Investigate the plasma velocities of 
laser-tarcet interaction using electrostatic probes tut with 
different target materials (metals) . An increase in plasma 
ion mass shculd markedly reduce their velocities. If this 
is pursued fcr research, the laser-target interactions of 
the cxides cf the target metals should also be investigated 
to resolve the guestion of how much the oxides contribute to 
the plasma picbe signal. 

(2) A concurrent project for further research is 
tc use alternate plasma diagnostic technigues for 
examination cf the laser produced plasma characteristics. 
Laser-plasma research at the Naval Postgraduate School has 
only beet dene using electrostatic probes, magnetic probes, 
pressure probes and X-ray photodiode prebes. Curved 
capacitor erergy analyzer probes, spectrcsccpic techniques 
and using a second laser shining through the plasma all 
offer very irteresting alternatives for collecting the laser 
produced plasma characteristics data £45]. Huddlestcne and 
lecnard [25] and Swift [53] shculd be consulted as reference 
tc hegin research in this area. 

(3) An investigation of the comments made by 
Porteus [41] and Isenor [27] about the relationship of 
background pressure to plasma ion velocity would make a good 
experimental series for thesis research. 

(4) Research the fluid mechanics of the laser 
induced tarcet crater formation. This is an area fcr some 
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It would ccntrikcte some 



gccd, tbcuctt prcvcking study, 
necessary atd useful knowledge about the laser target 
interacticD . 

(5) Another project for someone who is inteiested 
in boundary value problems and is good at computer analysis 
would be to take one of the simple, analytically unsclvable 
plasma models such as the NRL models in reference [37], 
£38], and [3S] and set it up for computer analysis. A good 
course to take if this project is undertaken is MA 3243, 
Numerical Methods for Partial Differential Equations, a core 
course in the Meteorology program. Reference texts by Smith 
[5C], Ames [1], and Spiegel [51] would also assist research 
in this area. 

Another project, more of a systems nature, 
would be to investigate the systems aspect of having 
available a reasonable annual budget (350K) for plasma 
research. ihat direction should the spending proceed would 
be the hypothesis with an eye toward the Laser-Plasma 
Laboratory as a complete system including electrical power 
requirements, space availability, equipment advances, 
student population, and the total cost figures fcr the 
operation. If the Laser Plasma Laboratory at the Naval 
Postgraduate School is to continue to produce good results, 
the publicity from such a thesis would be a beneficial side 
effect. A good paper to start with for such a research 
project is reference [5]. A visit ro laboratories at NHC 
China Lake, Berkeley, Livermore, and Sandia would also give 
a student of systems analysis a look at what state of the 
art research equipment is available. 

(7) If the laser pulses are continued after a 
crater hole has passed all the way through the target, does 
the probe sicral change? Investigation of this question 
would give a more complete understanding of the laser-target 
interaction . 

(8) The statements made in this paper about the 
energy velocity groups possibly being associated with the 
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aluiinum icEizaticn should be investigated. If results of 
tlat sort cat he derived using target material other than 
alurinum, tie hypothesis would be supported by a more 
conclusive argument. 

Alsc^ the relationship of the velocity ratios needs 
further wcrk. An explanation for the apparent gecnetric 
symmetry of the ratios was not discussed here. Further work 
with this theory should try to find a descriptive model to 
resolve it. 

(S) The plasma velocity as a function of mass was 
attempted, hut, due to time limitations, was not completed. 
There appear to be signals grouped at various velocities in 
Figure 1i and another possible explanation for the existence 
of the various velocity groups, besides the ionization 
theory, could be a mass theory. The slow signals may be 
produced hy heavy atoms (AL,0) ard the fast signals hy light 
particles (F.) , The suggestions by Isenor [27] and Pcrteus 
[41] that chcmher contaninents cause the large early-time 
peaks is further impetus to a theory dependent on particle 
masses. 
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APPENDIX A 



THE LASEE SYSTEM 



The Naval Postgraduate School Laser-Plasma 
Nd:glass laser system xas used for the experiment 
in this thesis. The system is composed of a Nd;gl 
a vacuum chamber, a system for plasma detection 
necessary ancillary eguipment to support each subs 
this thesis the floating double electrostatic prob 
was used for plasma diagnostics. A complete expl 
the flcatinc double probe used is contained in App 

A. THE lASEE SYSTEM DESCRIPTION 

A blcch diagram of the system components us 
thesis is pictured in Figure 17. The Nd:gl 

subsystem is centered around the KCHAD K-1500 
necdymiuii ir glass laser [57]. The laser output 
controlled by varying the voltage applied to the 
(K-1) and the amplifier (K-2) . The energy rang 
laser beam is 2.5 to 15 J. and the power of t 
approximately 540 MW. The pulse half width val 
nsec. The laser pulse shape is pictured in Figure 
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FIGURE 17. BLOCK DIAGRAM OF THE LASER SYSTEM COMPONENTS 
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FIGURE 18. THE LASER PULSE SHAPE 
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a 5C ca focal length lens was used to focus the laser 
keani dcwn tc a C.01 cm diameter, circular spot si 2 e. The 
flux density delivered under these conditions is akcut 

10 2 

5/ca . A, general description cf the KORAB Nd:glass 
laser physics is given by Maiman [34], The system installed 
at Naval icstgraduate School Laser-Plasma Laboratcry is 
detailed ky Eavis [16], 

tack laser pulse passes through the beam splitter 
which presides signals for two functions. The first purpose 
the beam silitter serves is tc reflect an electrcmacnetic 
signal tc the magnesium oxide photodiode detector (?-Cl) . 
That signal is fed to a Tektronix 564E storage oscillcsccpe 
which integrates it tc measure the laser pulse energy. The 
ether function of the beam splitter is to provide a trigger 
signal fer the Tektronix 7904 oscilloscope via a photodiode 
detector. That insures reproducible triggering of the 7904 
oscillcsccpe at the same instant of time. It also marks the 
arrival cf the flasma probe detection signal on every shot 
without losirg part of the laser pulse signal. Both the 
preke signal and the laser pulse can be displayed cn every 
laser shet. 

The laser pulse signal received by the MgO photodiode 
detector fren the keam splitter was much greater (about 10 
times) than the plasma signal from the probe. To scale it 
down tc a signal that can be displayed with the probe signal 
in the sane trace, a Tektronix bedel 3752 attenuator kcx was 
used with a 5C ohm resistor at either end and connected on 
one side tc the Tektronix 7704 oscilloscope and the other 
side tc the coaxial cable from the MgO photodiode detector. 
It acted as a voltage divider for the terminated signal, 
reducing its value (see Figure 19.) There are 10 
attenuators cf various resistances for ready use. 
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FIGURE 19. THE SIGNAL ATTENUATOR 
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If none fits the experiment r eguir ements , it is a simple 
matter tc change the resistance in the attenuate! by 
relieving the screws in the cover and soldering in a 
different resistance. 

Whet the laser beam hits the target in the vacuum 

7 2 

chamber with a flux density greater than 10 W/cm , a plasma 
plume is fermed. This plume was the subject under 
investigaticr in this research. 

The laser rcc cooling system cools the rods in the K-2 
amplifier ard the K-1 oscillator. There is a small 
refrigerator, a water pump, a water filter, two 
thermometers, a thermostatic temperature control and a 
reservoir cf deicnized water in the cooling system. 

There are four power devices associated with ccntrol 
of the laser beam. A 10 kilovolt power supply provides 
power tc the K-2 amplifier. A 5 kilovolt power supply 
provides power tc the K-1 oscillator. The shutter ccntrol 
is asscciatec with the Pockels cell ccntrcl of the laser 
pulse shape and the small 2000 volt HV-1565 power supply 
energizes the KgO photodiode detector. 

A CW EeNe laser provides visual alignmenx and target 
pcsiticnirg reference for the laser system. The laser 
system is diagrammed in Figure 20. 
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FIGURE 20. THE Nd:GLASS LASER SUBSYSTEM 
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Fct th€ detector subsystem associated with the laser 
(see figure 21) there are various diagnostic techrigues 
available for investigating the laser produced plasma. In 
the experime itatio E relevant to this thesis, a floating 
double probe (see Figure 22) was used for plasma deteotion. 
Ibe laboratory is also equipped to do researoh with macnetic 
probes, single filament eleotrostat ic probes and X-ray 
detectioE devices. The two sizes of plates for the top port 
in the vacuum chamber are tapped for 0,5 cm and 0.25 inch 
diameter probe shafts. Other detectors can be designed to 
the system. 

The plasma detector subsystem is made up of the 
photodiode to detect the initial laser pulse, a phctcdicde 
to trigger tfce oscilloscope on which the plasma energy 
signal appears, a Tektronix 5648 storage oscilloscope to 
measure the laser pulse energy, a Tektronix 7704 
cscillcsccpe to display the laser pulse shape and the plasma 
detector sicral, the plasma detector (an electrostatic 
floating double prcbe) , the detector power supply and 
coaxial cable delay lines for time signal coordination. The 
delay line from the probe to the 7704 oscilloscope allows 
for signal tine lac by a fixed amount to match the relative 
time dependence of the laser pulse, the trigger pulse, and 
the electrostatic probe signal. The signal speed through 
the coaxial cable is two-thirds the speed of light. A 
schematic of this subsystem is shown in Figure 21- 
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THE DETECTOR SUBSYSTEM 
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FIGURE 22. A FLOATING DOUBLE PROBE 
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ah€r the lasec beam is fired, the firing kej also 
energizes the camera on the oscillcscof a. A Tektronix flcdel 
3 automatic shutter control and camera mount and a Eclctcid 
camera {Eclcrcid 410 film) were used for the picture taking. 

A straj triggering signal will cccasicnally arrive at 
the Qscillcsccpe, triggering the camera early. Nothing 
appears cn the scope trace hut a straight line. If this 
happens tco cften, check the oscilloscope sensitivity and 
slope settings. They were a very fine, easy to disturb 
adjustment during these experiments. The time constant for 
any cable can he determined on the Tektronix 519 
cscillcsccpe hy the following procedure. 
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JijES of Coaxial Cables 



Plug in the scope and turn it "ON". 
Connect the HATE GENESATOB +HAIE 50 
into the INPCIS SIGNAL 125 Q terminal through a 
connecter . 



f 

Q terminal 
T50/N125 2i 



On the THIGGSH panel, set the lEIGGEH SCU 
HATE GEK. The room should be in ssmidarkness to re 
purple oscilloscope trace. 

Set the HATE GSKEBATCH panel MUITIPLIEH 
tc X100 and the CICLES/S3C X «13LT. switch between 6 an 
On the TfilGGEH panel, set the GAIN tc N 
the FCNCTION switch tc PULSE and the PULSE AHPLITUDE 0 
switch to the black arrow indicaticn. 

Set the TIME BASE panel NANOSEC/CM switch 



and the 


EELA Y knob to 


the appr 


cxiaate center 


cf 


scale . 




Insert a 


coaxial 


cable connect 


or 


tee a 


INPUTS 


SIGNAL ^2SQ 


line 


on the end 


cf 


the 15 



BCE to 
ad the 



switch 
d 12. 
OHMAL, 
F SYNC 

tc 10 

t the 
C/N125 



adapter. 



Beccnnect the +BATE 50 Q lead to one side of 
the tee and the cable to be measured tc the other side of 
the tee. 





look on 


the scope 


leading edge 


to the 


signal is 


to go up and 


back on 


the cable 


is the cable 


's time 


delay. 



trace . 


The 


time 


from the 


the 


time 


it takes for a signal 


so 


divide it 


by two 


and that 
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C. IHE ViCOCe CHAeSER 



Ih€ vacuum chamter subsystem is the final part cf the 
laser sjsten (see Figure 23) . It consists of the vacuum 
chamber, a gas bleeder system, the target and target hclder, 
an cil diffusion pump, a liquid nitregen cold trap, a 
mechanical fere puop, a vacuum ctamber thermocouple gauge, 
an icn gauge in tie diffusion pump, a thermccouple gauge in 



the fere 


fuirp vacuum 


line, an 


ion 


gauge in the 


vacuu m 


chamber 


and a 


mete 


t to read 


the 


pressure at each 


gauge 


(GMA-201 




and 


IONIZATION 


GAUGE CCSTBOLLEE) . 


The 


f cr epump 
- 13 


brings 


the 

-3 


system vacuum 


pressure down tc 


about 



10 torr. At 10 torr the cil diffusion pump takes ever 
and can reduce tie pressure in the vacuum chamber furtier tc 



, 7 

arcut 1C terr, Ihe liquid nitrogen trap areve rhe 

diffusicr pump acts as a condenser for rhe hot oil drcplets 
coming up frem the heater. Ihe diffusion pump exterior is 
ceded bj an open ended tap water cooling system that drains 
off into a well outside of the building. 

The vacuum chamber can be isclated from the vacuum 
pumps by the chamber gate valve. The vacuum chamber has 
five large side ports for instrumentation and observaticn, 
twe small sice ports, and another large port on top, also 
for diagccstic instrumentation. The target (usually a metal 
disc) is suspended in the vacuum chamber at the end cf the 
threaded metal red wirh a knurled handle outside cf the 
chamber. The rod passes through the vacuum chamber wall in 
a vacuum fitting to facilitate the experiment operacion by 
expesing a clean target surface tc the laser team 
periodically without breaking the chamber vacuum. 
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FOREPUMP 




FIGURE 



23 



THE VACUUM CHAMBER SUBSYSTEM 




DIFFUSION 

PUMP 



FIGURE 



24 . 



ROUGHING VALVE OPERATION 
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Ih€ vaiicus gauges (twc thermccouple and twc ion 
gauges) ncnitcr the system fressure at various points in the 
system (see figure 23) . Each gauge is attached to an 
electronic neter to read the pressures at the mcnitcring 
points . 

The fcrepump was used to rough pump both the vacuum 
chamber and the diffusion pump through a system of valves 
under the vacuum chamber (see Figure 24). Cnee the 
diffusion pump has been cut in tc the vacuum chamber the red 
roughing valve must remain closed to prevent hack pressure. 
The procedures in the Laser light Off section cover this 
idea thoroughly. 

A recent addition to the vacuum chamber subsystem for 
fire safety and prctecticn of the diffusion pump are two 
relay switches. A thermal switch was added tc prevent the 
cil diffusicE pump frem ever heating and burning cut the 
heater ceil, and a flow switch was added tc shut the oil 
diffusicE ptmp off if the flew cf the coding water to the 
pump is icterr uptad . The electric circuit schematic fer the 
switching devices is shown in Figure 25. 

The laser system safety aspects must be well 
understood hj the researcher. Tc remember the standard Navy 
safety precautions for working around high voltages is a 
definite way tc avoid an accident. familiarity and 
overlooking common sense for the sake cf expediency is just 
as cangercus around a high powered laser system as it is 
aboard ship. All cf the safety interlocks in the world will 
net protect a technician from improperly grounded or shorted 
system safety features. A very good summary cf the safety 
aspects cf lasers is contained in Chapter 7 cf Eeady [42] 
and should he read by all laser system operators often 
enough fer their lessens to become second nature. 

The laser system operating proced">res are included in 
the fcllcwirc sections to give future students ready access 
tc lessens already learned about the eguipment operation. 
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lASEE OEiEATING FBOCEDURES 



1* Nc^Glass Laser Li ght-Off Pr ocedure s 



1. Ecrefunj 

Close siai shaped vacuum vent. 

Push fcrepump "SIART” button. 

Start the air conditioners to dehumidify the rccm by 
switching 411c and #9b in the circuit breaker panel. 

2. Ion Gauges 

lurn cr the power switch to GMA-201. This is a 
thermocouple gauge to measure the chamber pressure. Select 
station 41 and this meter will operate in an automatic mode. 

Turn on the power to 8G-3A. The thermocouple gauge on 
this meter is set at T.C. #1. It measures the forepuiiip 
rctcb vacuum. 

The degas switch sends a current through the ion tube 
tc expel ahscrbed gas from the electrodes. The system has 

— c 

tc be at 1C torr cr less to use this feature. 

The current set, current adjust, and read current 
indicators let you observe the ion gauge grid parameters 
frcir time tc time. 

3. Diffusicr Pump 

Turn on the cooling water. 

Don’t turn on rhe diffusion pump until the forepuap has 
been operating for about twenty minutes and the pressure in 
the diffusion pump system is less than one hundred microns 
because the diffusion pump will not pump dcwn in greater 
pressures. Also operating at excesive pressure destroys the 
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fuap. Tfc€ diffusicn cil overheats and evaporates anc 
heating ceil turns out. 

Make sure the chamber gate valve is closed. 

The Hack valve under the chamber should be open. 

The red valve under the chamber should be shut. 

Push the diffusion pump ''START" button. The red "ON" 
light on the metal relay box below the button will come on. 
The relay hex is wired to a temperature sensitive thermal 
switch on the diffusicn pump and a pressure switch on the 
cooling water line. The thermal switch will shut off the 
diffusion pump if it overheats and the pressure switch will 
shut eff the diffusion pump if the cooling water line loses 
tack pressure, ie. if the cooling water is shut off before 
the pump cr the coding line starts to leak. 

4. Ion Gauges 

When the diffusion pump has been on for ten minutes, 
the pressure in the diffusicn pump system will be 
sufficiently low to turn on the ion gauge in the diffusion 
pump system without burning out the filament. A vacuum of 

-3 

abcut 1C terr or less is required cr an automatic meter 
feature' will shut it off. 

-4 

Set the pressure multiplier selector tc 10 

Adjust the zero cn the ionization gauge meter. 

Push the filament "ON" button. The red light should 
ceme on and the ion gauge tube filament at the diffusion 
pump should start tc glow. 



5. Diffusicn Pump Control 
Let the diffusion 
minutes. while recording 
ferepump thermocouple pres 
eff the electron tube ion 
glass in the system tend 



pump operate fer about 
the ion gauge pressure an 
sure, note that whenever ycu 
gauge detectors, the wire 
tc heat up. There is a tem 



twenty 
c the 
light 
s and 
porary 
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rise in fiessire dee to the degassing at the solid surfaces 
in the sjsten «hich implies increased pressure. 

When the pressure in the diffusion pump chamher has 
-s 

reached 10 torr (in about twenty minutes) , add the liquid 
nitrogen to tlie cold trap. The pressure will drop two more 

-7 

orders cf magnitude to about 4x10 terr. This step is 

-5 

optional if 1C is an acceptable chamber pressure. 

6. Opening the Chamber to Vacuum 

- 4 

Turn the ion gauge scale to 10 , because when the gate 

valve is opened there will be a rush cf pressure tc the 
diffusion pump which will cause the detector tc impulse to 

-4 

about 5x1C terr. 

Make sure the chamber vacuum is tight. 

Close tbe black valve belcw the chamber. 

Crack the red valve when satisfied that there are no 
leeks. 

If the pressure, as indicated by the thermcccuple 
gauge, is decreasing, then the chamber is probably vacuum 
tight. If a very porous or very effusive material (tape, 
paper, fresh paint, etc.) is put in the chamber, it may take 
a very Iccg -time tc get a reading below cne cr two hundred 
microns. 

When the pressure in the chamber reaches sevent y-^f ive 
micrens cr less, close the red valve. 

Open the black valve. 

Crack tie chamber gate valve epen. 

Check tc see that the pressure in the chamber is 
remaining at cne hundred micrens cr less and gradually open 
tbe gate valve. If the pressure does not decrease, close 
the gate valve and return tc the beginning cf this 
procedure . 

Continue to record the ion gauge pressure, ferepump 
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pressure, anc chamter pressure for about ten minutes to 
insure tfce sjsteir is vacuum tight, 

7. Chamter let Gauge 

This detector is easily burnt out and bard to replace 
sc it is imperative that the system pump fer about fifteen 
minutes before turning it on and that it be allowed tc cool 
off for fifteen minutes before the chamber is returned to 
atmospheric pressure. Diffusion pump pressure should be of 

- e 

the order of 5s10 torr or less when the chamber icn gauge 
is turned cn. 

Set the ionization gauge controller power switch to 

«CK''. 

Wait abcut thirty seconds for it tc warm up. 

Select pressure on the function switch. 

Push meter zero-ro-check and adjust as required. 

Select emission current 3.2-10 mA. 

New select "ADTO'' on the power switch. This 
automatically lets the meter circuit select the vacuum 
scale . 

Push tie filament ••ON” button. 

The initial reading of the gauge meter will be ir the 
-s 

range of 1X1C terr due to out-gassing. 

The final pressure should be reached in five minutes at 
■. « 

about 5x1C torr. The limit of this vacuum system is a 

-6 

chamber pressure of 1.5x10 torr and diffusion pump 

..8 

pressure at about 9x10 torr. 

The chamber is now in vacuum and ready for the laser 
experiments. Note that the laser shots create plasmas which 
increase the pressure readings in the system briefly. The 
chamber icn gauge can automatically compensate for tho 
change by shifting scales but the diffusion pump icn gauge 
does not have this feature sc it is good procedure tc set it 
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tc the 1C scale when shooting a target in the chamfcer. 



8. Ihe laser 

First turn on the continuous He-Ne laser to see 
the Nd laser beam will be going. Turn the switch 
bach of the Ee-Ne laser to *'0N". Sait three seconds 
push the switch to "START". Be careful net to no 

laser bj rcugh handling when switching because 
cpticallj alligned. 



where 
at the 
then 
ve the 
it is 



9. The laser Cooling System 



Turn cn the right switch 


at the 


w? 


ter pump. 


The 


coding water is special 


deionized 


water. 


DC 


not use 


tap 


water in this system. 


Check tc 


see that 


the 


water bubbles 


through the transparent 


tee over 


the re 


ser 


voir. If 


it 



doesn't, there are problems in the cooling waver systen. 

Turn on the left switch to start the refrigerator. 
There is an automatic temperature ccntrcl just under the 
transparent tee that can be adjusted with a screwdriver to 
ccntrcl the ceding water temperature. 



10. Laser Eewer Supplies 

Eush the shutter control power butter. It will light 
up when pewer ccaes cn. 







Switch 


cn 


the five 


kilovolt 


?o 


wer 


supply 


• 


The 


"B 


ONE" 


lig 


ht 


wi 11 c 


CHS 


cn when po 


wer is 0 


n. 


















Switch 


cn 


the ten 


kilovolt 


po 


wer 


supply 


• 


The 


"D 


DMP" 


lig 


ht 


will c 


CHS 


on when po 


wer is 


on . 


Let po 


wer 


su 




lies 


war 


It 


up f cr 


live 


minutes b 


efere us 


ing 


the 
















Switch 


"AC 


ON" switch 


tc the 


HV- 


1565 


power 


su 


pply 


• 








Switch 


"HV 


CN" switch 


to the 


HV- 


1565 


pewer 


su 


pply 


• 








Sait fc 


r ten minutes 


then che 


ck 


the 


temp 
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at the K-i anplifier. Safe laser cooling water tempe 
is about 20° C. Kcnitor the temperatures cccasionall 
the cooling water is inadvertently overlooked and not 
CD, the temperature will gc up on each shot. I 
happens CC NCl IQfiN ON THE CCCLING WaTER, shut the 
system cff for at least two hours so the red in the 
will cool off naturally. If the cooling water is turn 
the sudden heat flow will break the glass red, mini 
laser (this has happened) . 



11. Experimental Eguipment 

Turn on all cscillcscopes, cameras, expe 
egeipment pewer supplies, and all other eguipment f 
experiment tc be run. 



12. Laser hiring 

Blow cff all optical surfaces of the laser with a 
bull befere starting. Dust and dirt will tend tc 
pitting cn the optical surfaces whenever the laser is 
Cc not blew cn any optical surface with your b 
Heisture will deposit and, like dust, will tend to pi 
surface . 

There are two methods of firing the Nd laser, Icc 
remete. Ihe local method is from the 10KV power supp 
the remete is from a specially constructed control bex 

Hhec ready to charge the power supplies to she 
laser, anncunce the fact to all in the rccn as a 
precauticr . 

Eut your safety glasses in place. Do net proceed 
every one has cn safety glasses. 

Hher everyone is ready, push the ’'CHARGE" button, 
vcltage will start coming up on the meters at the 
supplies. Chiming bells will toll at the rate of abcu 
per secerd urtil the system is completely charged. Hh 
power supplies are ready to fire the laser, a cent 
mcrctone whistle will sound. The ready light will c 



nature 
j. If 
turned 
f that 
whole 
laser 
ed on 
ng the 



riment 
or the 



r air 
cause 
fired, 
reath . 
t the 

al and 
ly and 

• 

ct the 
safety 

until 

The 
power 
t two 
en the 
inuous 
cme on 



101 



-ft! in 
' tk a i 

- 

-1 *1 iif 

tu« 

f> £ ^ tf 

' z » 1 •£ 

1 . -f f 

)l# 



..I 



: ilOz 

( 4i|ti| 
'Oil all 
ziilalofi 

A_ 

• arlxnjz 

^ l 

t*^3 atf# 



\3frft 

T . 



♦ 



CD tke r€BCt€ control and the "FIP.£‘' light will come cn at 
the 10KV and 5KV pcwer supplies. 

If fcr sene reason the shot is to be aborted, push the 
"CDUE” button cn the 10KV power supply or the remote control 
box. Do not ”EOMP'‘ or ’’FIRE" before the power supplies are 
fully charged because it causes excessive currents in power 
supply control relays and tends to burn them out. 

To shoot the laser, push the ”FIEE” button on the 10KV 
pewer supply or on the remote control box. 

The time between consecutive shots should be sixty 
seconds or mere to allow time for the glass rod in the laser 
to recover a constant temperature gradient. 

Check the temperature over the lOKV power supply air 
vent cccasictally to detect any overheating. 



2* Shut Eewn Proce dure s 

1. Experimert Eguipment 

Shut off all oscilloscopes, cameras, experiment power 
supplies, and other experiment equipment power. 

2. The Laser 

Turn off the 10K7 power supply. 

Turn off the 5KV power supply. 

Turn off the shutter control power by pushing the power 
button. The light turns off when power goes cut. 

Turn off the HV-1565 high voltage. 

Turn off the HV-1565 AC. 

Turn off the continuous He-Ne laser. 

Let laser cooling water flow for ten minutes after the 
last laser shot, then turn off the two laser cooling water 
control switches. 
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. Vacuua Cbanber 

Ieav€ GKi-201 on indefinitely. It will not damag 
E€t€I. 

Turn off the chamber ion gauge filament. 

Turn off tte chamber ion gauge filament power. 

Close tie chamber gate valve. 

Let the ion gauge filament cool down fcr about £ 
minutes tefcre venting the chamber. Failure to let it 
tends to ttrc cut the ion gauge filament when the cham 
vented. Ihe vent is located above the red valve unde 
chamber. Vent the chamber if desired. 

4. Eiffcsicr Eump , 

.4 

Select 1C cn the pressure multiplier switch tc 

Turn off filament to the diffusion pump ion 
power can le left cn continually to EG-3A without da 
the meter. 

Eush ”STCE" button on the diffusion pump. 

let the diffusion pump cooling water continue to 
fcr about twenty minutes. This is an open end c 
system that dumps cut in the well outside of the bu 
inside of the theta-pinch power plant fence. The dif 
pump and the cooling water can be left cn indefinite 
desired while research is in progress. 

5. Wrap up 

Turn off the air conditioners at circuit b 
switches 111c and #Sb. 

Leave tte forepump on continuously. To shut 
when desired, push the "STOF" button and open the 
shaped vent valve. 
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aPF£NDIX B 



THE ELECTBCSTATIC EBOBE 



used to 
plasma p 
voltage 
plasma, 
the elec 
isclatic 
characte 
characte 
shi eld 
Debye s 
pctentia 
lengths , 



Assuming 

indicate 

- 3 

[cm ]. 



Ihe electrostatic floating double probe signal was 
determine the time at which the laser produced 



ulse arri 


ved 


at the 


probe location. 


When a 


DC 


is applied 


between 


two electrodes 


in a neu 


tral 


the ions 


are 


attracted 


by the negative 


electrode 


and 


trcrs by 


the 


positive 


one. The charge 


migration 


and 


n crtates 


space charge 


electric fields. 


changing 


the 


r of 


the 


plasma 


neutrality. A 


f undaae 


ntal 



ristic of the behavior of a plasma is its ability to 

itself from the electric fields in and around it by 

hielding. The sheath thus formed absorbs the 

1 difference exponentially within a few Debye 

X , where, 

D 

1/2 

X = 69(T /n ) [meters], 

C e € 

T = [OR], 
e 

-3 

n = [m ]. 
e 

n and n are approximately equal. Table 7 
i e 

11 13 

s that n varies from about 1.08X10 to 7.3MX10 
e 

T , oaloulated for this experiment by electrostatic 
e 



104 



f 



smaif 

, ca»s£ 

- f a 



•i»ld 

f 

.1 totfsl* 



•leil. 



.fit 



♦ 



picb€ theory, was about 50 eV. On the same experinental 
apparatus, Slewchuk [49], using X-ray diagnostics, get an 
electron temperature at the target surface of 162 eV. 
Williamson £61] also used a temperature of 162 eV in his 
experiments using magnetic probes at the target surface. 
Biceks [1C] in earlier research used an average temperature 
cf 100 eV fer his calculations. A plasma temperature of 50 
eV is consistent with Shewchuk's value two cm from the 
target surface. 

To use the electrostatic probe theory in discussions 
of the plasma characteristics, the probe radius must be much 
greater than the Debye length which is the case for a Debye 

- 3 

length cf about 10 cm at two cm from the target surface. 

Using n = 7.34X10 cm and T =50 eV, the minimum Debye 

e B 

-♦ 

length at twe cm from the target was 1.94X10 cm. Using 

11 -3 

n = 1.1X1C cm and T = 162 eV, the maximum value for 
e e 

the Debye length in the plasma right at the target surface 
was 2.85X1C cm. Then the Debye length is cf the erder of; 



-2 

2.65X10 > X > 

D 



- 4 

1.94X10 cm. 



-2 



Ihe probe radius was 2.54X10 c 

egual to the greatest possibl 
i.e. at the target surface, 
approximately 130 times larger t 
measured at two centimeters from 
An advantage of the double 
types cf probes is that the maxi 



m w 

e va 
Th 
ban 
the 
pro 
mum 



hich 


was 


apprexi 


lue 


of the 


Debye 1 


e 


prebe 


radius 


the 


miniou 


m Debye 



target surface, 
be system over the 
current drawn is 



irately 

ength, 

was 

length 

other 

always 
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liirited by the ion densities because the cuirent is limited 
with the ion saturation current. In the double probe theory 



developed 
[amp], is 



where : 



for a static plasma the saturation ion current. 



related to the density by the relationship 

1/2 

n = 1.671 (m /kT ) /eA 
i + i e p 

2 



A = surface area of the probe [cca ] 
P 

- 3 

n = icn concentration [cm ] 
i 

m = ion mass [gm] 
i 

e = particle charge [ccul] 

kT = electron temperature [erg] 
e 

and 1.67 is a unitless ccnstart. 



I 

+ 



The electron temperature can be deduced from the slope 
of the I versus 1 probe characteristic curve at V = C (see 
Figure 5) , using the relationship: 

.. 1 

T = (el /2k) (dl/dV) [°K] 
e + 

where I is the saturation ion current in amps. 

+ 

For the case of a streaming plasma the ion current is 
given [ 1 0 [ by : 

I = n ZeA V 
+ i p LP 

where: n Ze = average ion charge per unit volume 

i 

1 = plasma stream velocity [cm/sec]. 

LP 

Or acre accurately; 

I = n Ze(1 f s)A V 
+ i p LP 

where the coefficient, s, is included to correct for the 
secondary electron emission from the probe. The values cf s 
range from 0.0 to 0.5. For tungsten probe filaments s is 
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apprcximatelj egual to 0.5. 

For the experiment performed, the average ion 

densities in the plasma plume, when the probe was plaoed 2 

cm from the target, were calculated and arranged in Table 9. 

I = V/B = n Ze (1 + s) A V 
+ i p LP 

where: V = probe response voltage [volts] 

A = 2 jr rh, the exposed probe surface area 
P 

h = 0.32 cm, the whisker length 
-2 

r = 2.54X10 cm, the whisker radius 

B “ total circuit resistance in ohms 

n = average ion density at 2 cm 
i 

solving for o using, 
i 

A 

P 
B 
s 

Ze 

then n 

i 



- 2 2 

- 5.1X10 cm 

= 1 ohm 
= 0.5 

.19 

= 1.5(1.6X10 ) coul 

19 -3 

= 5.45X10 V/V [cm ]. 

LP 
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FIGURE 26. 



DOUBLE ELECTROSTATIC PROBE CIRCUIT DIAGRAM 
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lable 9. Average Icn Densities (r ) 

i 

. 3 



1 fcir/sec] 


V [vclts] 


n [cm 3 




IE 




i 




6 




1 3 




4. C0X10 


3.0 


4.09X10 




6 




1 3 




5.57X1C 


7.5 


7.34X10 




7 




1 3 




1,52X10 


2.8 


1. 00X10 




7 




1 2 




3.71X10 


0.7 


1.03X10 




8 




1 1 




1.01X10 


0.2 


1.08X10 




Electrostatic prcfaes 


are actually quite simple 


to 


ccEstruct (see figure 22) . All 


that is needed are 


two 


suitatle fcire whiskers, some method to insert them intc 


the 


vacuum cbamter, and an 


electronic 


circuit to display 


the 


signal or an cscil Icscope. 


The probe 


circuit used in 


this 


experimentation, designed 


by Brooks 


[10] is shown in Figure 


26. Additicnal electrostatic probe construction techniques 


are covered in Chen [25] 


pa.ges 194 


thru 199. Table 10 


can 


be used to check the probe 


circuit 


integrity. Refer 


to 



Figure 2€ fci the test points. 
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lable 10 



Electrostatic Prote Circuit Verification lafcle 



Resistance Actual Probe 
Reading Resistance 

Ccnnecticn Good Readings 

(ohms) (ohms) 

A to I Q 0 

E to J 0 0 



E 


to 


E 


>100K 


C 


to 


E 


> 100K 


E 


to 


E 


0 


G 


to 


E 


>5K 


E 


to 


C 


1 


E 


to 


B 


0 


E 


to 


C 


1 


Problem 


ai€as 


to watch f 


Ity are 


• 







IM 

2H 

0 

10K 

1 

0 

1 

r when the jrcbe appears to 



(1) The BNC connectors at HG and DC are no 
well irsulated from the probe chassis. If they 
cut, they should be replaced with BNC connector 
have at insulating plastic ferrule fixed to 
connecter tc insulate it freza the chassis. 

(2) The 4 microfarad capacitor went out 
this experiment. It is not a standard size sc 
were stocked in the laboratory in case another fa 

(3) The probe whiskers AI and EJ can tre 
very easily and great care must be exercised 
hanclinc the probe. Only insert and remove the 
from the vacuum chamber when necessary. Handle 
little as possible and the probability of breakag 
be Diciaized. 

(4) After forty or fifty laser shots the 
whiskers get covered with molten aluminum glo 
They should be carefully cleaned off with a knife 



t very 
ground 
s that 
the 

curing 

extras 

ils. 

ok off 

when 
probe 
it as 
€ will 

probe 

hules. 

blade 
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k 




i 
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at ccnverient intervals during the experimentation. 



& C . 2 
outside cf E 
The over a 
tte chassis, 
the shaft 
Very Icng pr 
inserting a 
hclder acc a 
support for 
using chentis 
structural t 



5 inch diameter metal tube slipped ever the 
recks [10] probe will provide added strength. 
31 length of the probe, from the whiskers up to 
should not exceed 20 cm in length. The shorter 
the easier it is to handle it without breaking, 
ches have a tendency tc break easily when 
nc extracting them from the vacuum chamber probe 
re hard tc store and work around. Extra shaft 
existing probes with glass shafts can he had by 
try laboratory stands, clamps and solid metal 
ubes . 
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APPEHDIX C 



IHB COHFUTATIONAL ANALYSIS C5 EATA 



Ic f^cilit 
Cfatrch CCDfUt€r Ce 
Th€ follcwicg re^€ 
ptcgranis: Gcttf 
oscilloscope phct 
analy2€d and trans 
cata taker ficw ea 
consisted cf the f 



ate experimental data red 
nter IBM 360 ccmputer sys 
rences were used to wri 
ried [22] and Shaap 
cgraph data from each 
ferred to data cards (see 
ch oscilloscope prone s 
cllcwing list of variable 



uction the 
tern was emp 
te the CO 
ine [ ^6 ] • 
laser she 
Figure 27) 
ignal phct 



a, H. 
Icyed. 
aputer 
The 
t was 
, The 
cgraph 



(1) The picture number, called matrix DESIG, given 
in a double precision form is XXXX.YY22AA, where XXXX 
is the local time of the picture, YY is the day cf the 
month that the picture was taken, ZZ is the menth, and 
AA is the last two digits cf the year^ 

(2) The laser shot number (the number one being a 
fresh target surface) was called matrix SHOT fer data 
input and later in the program SHC for additional 
sorting capability. 

(3) The distance of the detecticr probe from the 
target surface was called matrix D for input and . DI3 
fer sorting. 

(4) The probe DC voltage bias was called BIAS. 

(5) The time of the peaks in the oscillcsccpe 
pictures was called matrices T1, T2, and T3 for data 
input because some photographs had as many as three 
readable peak values. Time was changed to matrix T for 





Id) 

&i sld> 

stzaaa 



'f 

i)iiilli 






scrticg acd ccmputing operations. 

(6) The natrices PK1, fK2, and PK3 were used to 
input the peak values of the electrostatic probe 
signals. The matrix PK was used later in the pregram 
fer ccmputing. 

(7) The laser pulse energy was designated matrix E 
for input and changed to EN fer ccmputing. 

(8) The measurement errors fer time and peak probe 
voltage values were recorded as matrix ERT and ESVOL 
respectively. These were to allow the experimental 
error bars to be computed and plotted. 

The typical data card was punched as pictured in 
Figure 71. h total of 376 pictures were included in the 
data deck fci machine analysis. Nearly as many pictures 
were discarded. They were unreadable because the time mark 
was net readable or because no peaks appeared cr the 
oscillcsccpe scales picked for the shots. 

The original pregram to execute data reduction was 
written in Icur sections. The first section was tc read in 
the data cards (see Figure 28) , the second tc locate all 
non-zero peaks and their associated time cf occurrence and 
change them tc matrices PK and T respectively (see figure 
2S) . The third section, also in Figure 29, operates cn the 
graphical peaks, ccmputing the distance the plasma had tc 
travel frem the laser crater tc the prete and using this 
distance to compute a velocity for each data peak. The 
remaining section (see Figure 30) was used tc graph the data 
points. A ccnputer library program DEAHP was used in this 
section tc produce the graphs included in the results of 
this thesis. The original program is listed here. 
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As scr 
sections ha 
Figure 31) i 
fcy shot nun 
served their 
stcrage spac 
Furthe 
data requir 
laser shct t 
Section six 
section cne 
to average 
sutroutire . 

A (1) = SI, a 
'standard dev 
ag to the pr 



e information was required from the data, 
c to te added to the program. Section £iv 
s a sorting routine to separate the data 
ter. Notice that the matrices that have a 
purpose can be used again to save compute 



e . 

r unanswered questions about the experi 
ed a subroutine to average the energy c 
c see if energy had any bearing on the re 
(see Figure 32) was inserted just fol 
and AVgNEB was a subroutine designed and 
the data pcint energies. AVENEB was iiste 
AVENEE takes matrix A with an initial 
verages the values A(1) to A (NB) and retur 
iaticn, SIG, and the average value of matr 
ccram. 
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FIGURE 28. COMPUTER PROGRAM SECTION ONE 
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N = C 

CC CCMINUE 
h = N-»i 

IF(CESIG(N).LE..01 ) GC TO 30 
If (T1(N).LE..01) GC TO 70 
FK(I)=FI<1<N)/BIAS(M 
T (I ) = T1 <^) 

^ - = 1 

GC TO 21 
CCMINUE 
1 = Ml 

7C CCNTINLE 

IFn2(N).LE..01J GC TO 80 
FK(I)=FK2(N)/BIAS(N) 

T (I )=T2 (^) 

N 2 = 2 

GC TO 21 
5C CCMINUE 
I = I41 

EC CCNTINLE 

If (T3{N ).LE..Oi) GC TO 90 
FK(I)=FK3(N)/BIAS(N) 

T (I ) = T3 (N) 

N 2 = 3 

GC TO 21 
EC CCMINUE 
1 = 1 + 1 
GC TO 9C 



SECTION 2. LOCATES ALL NON-ZERO DATA POINTS 



21 CCNTINLE 
SN(I)=£(N) 

SFC(I)=SHCT(N) 

C IS m = C(N ) + (. 01^ <SFQT(N )-!.)) 
IF(T( I) .LE..01 » GC TO 110 
VEL (I ) = (CISd )/T{ I ) )=» 1000000000. 
IF(N3.EC.3) GC TO 60 
IF(N3.EC.2) GO TO 50 
1 F(N3.EC.l ) GO TC AO 
2C CCNTINLE 



SECTION 3. COMPUTES PLASMA ION VELOCITIES 
FOR EACH DATA POINT 



FIGURE 29. COMPUTER PROGRAM SECTION TWO AND SECTION THREE 








k 



HE (2) = 1 
HE (3 ) = E 
HE (4) = 10 
I7E (10)=2 
HE (12J=1 
R 7 E ( 1)=.2 
F7E (2)=1.E7 

ECLIV^LENCE (TI T1 , FTB (5 ) ) 

RE/L^a TITK12) 

RE^C (5,300) TITl 
3CC FCFN)i7 (6/S8) 

C/LL CR/WP (NliTl,i 
ECLIV/LENC5 (TITl, 

RE/C (5,300 TITl 
C/LL CR/WP (N2,T2, 

. ECLIV/LENCE (TITl, I 
RE/C (5,300 TITl 
C/LL CP/WP (N3,T3, 

ECLIV/LENCE (TITl, 

PE/C (5,300 TITl 
C/LL DR/k<P (N4,PK1 
ECLIV/LENCE (TITl, 

RE/C (5,300 TITl 
C/LL CR/WP (N5,PK2 
ECLIV/LENCE (TITl,f 
PE/C (5,300 TITl 
C/LL DP/WP (N6,PK3,V£L1 ,IT8,RTB ) 
lie CCN7INUE 
S7CF 
ENC 



.6I/S,ITB,RTB) 
rRTB(5) ) 

,SFGT,ITB,RTB ) 
,RTB(5 ) ) 

,E,ITB,RTB) 

) R T B ( 5 I ) 

.,C,ITE,RTEJ 
,PT8(5 ) ) 

i,St-Gl,ITB,RTE ) 
rRTB(5 ) ) 



FIGURE 30. COMPUTER PROGRAM SECTION FOUR 
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FIGURE 31. COMPUTER PROGRAM SECTION FIVE 
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FIGURE 32. COMPUTER PROGRAM SECTION SIX AND SUBROUTINE AVENER 
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